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FOREWORD 

This  document  is  Volume  IX  of  the  Interim  Report  series  for  the  Passive 
Nosetip  Technology  (PANT)  program.  A summary  of  the  documents  in  this  series 
prepared  to  date  is  as  follows: 

Volume  I - Program  Overview  (U) 

Volume  II  - Environment  and  Material  Response  Procedures  for  Nosetip 
Design (U) 

Volume  III  - Surface  Roughness  Data 

Part  I - Experimental  Data 

Part  II  - Roughness  Augmented  Heating  Data  Correlation 
and  Analysis  (U) 

Part  III  - Boundary  Layer  Transition  Data  Correlation 
and  Analysis  (U) 

Volume  IV  - Heat  Transfer  and  Pressure  Distributions  on  Ablated 
Shapes 

Part  I - Experimental  Data 
Part  II  - Data  Correlation 

Volume  V - Definition  of  Shape  Change  Phenomenology  from  Low  Tempera- 
ture Ablator  Experiments 

Part  I - Experimental  Data,  Series  C (Preliminary  Test 
Series) 

Fart  IT  - Experimental  Data,  Series  D (Final  Test  Series) 

Part  III  - Shape  Change  Data  Correlation  and  Analysis 

Volume  VI  - Graphite  Ablation  Data  Correlation  and  Analysis  (U) 

Volume  VI T - Computer  User's  Manual,  Steady-State  Analysis  of  Ablating 
Nosetips  (SAAOT)  Program 

Volume  VIII  - Computer  User's  Manual,  Passive  Graphite  Ablating  Nosetip 
(PAGAN)  Program 

Volume  IX  - Unsteady  Flow  on  Ablated  Nosetip  Shapes  - PANT  Series  G 
Test  and  Analysis  Report 
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Volume  X - Summary  of  Experimental  and  Analytical  Results 

Volume  XI  - Analysis  and  Review  of  the  ABRES  Combustion  Test  Facility 
for  High  Pressure  Hyperthermal  Reentry  Nosetip  Systems 
Tests 

Volume  XII  - Nosetip  Transition  and  Shape  Change  Tests  in  the  AFFDL  50 
MW  RENT  Arc  - Data  Report 

Volume  XIII  - An  Experimental  Study  to  Evaluate  Heat  Transfer  Rates  to 
Scalloped  Surfaces  - Data  Report 

Volume  XIV  - An  Experimental  Study  to  Evaluate  the  Irregular  Nosetip 
Shape  Regime  - Data  Report 

Volume  XV  - Roughness  Induced  Transition  Experiments  - Data  Report 

This  report  was  prepared  Ly  Aerotherm  Division/Acurex  Corporation  under 
Contract  F04701-71-C-0027. '"Volumes  I through  IX  covered  PANT  activities  from 
April  1971  through  April  1973.  Volumes  X through  XV  represent  contract  efforts 
from  May  1973  to  December  1974.  Volume  X summarizes  the  respective  test  pro- 
grams and  describes  improvements  in  nosetip  analysis  capabilities.  Volume  XI 
presents  an  evaluation  of  the  ABRES  test  facility  in  terms  of  performing  ther- 
mostructural  and  reentry  flight  simulation  testing.  Volumes  XII  through  XV 
are  data  reports  which  summarize  the  experiments  performed  for  the  purpose  of 
defining  the  irregular  flight  regime.  The  analysis  of  these  data  are  presented 
in  Volume  X. 

This  work  was  administered  under  the  direction  of  the  Space  and  Missile 
Systems  Organisation  with  Lieutenant  A.  T.  Hopkins  and  Lieutenant  E.  G.  Taylor 
as  Project  Officers  with  Mr.  W.  Portenier  and  Dr.  R.  L.  Baker  of  the  Aerospace 
Corporation  serving  as  principal  technical  monitors.  Mr.  M.  J.  Abbett  was 
test  engineer  and  principal  Aerotherm  investigator  for  the  work  described  in 
this  volume. 

Thi.s  technical  report  has  been  reviewed  and  is  approved. 

E.  G.  Taylor,  Lt. , USAF 

Project  Officer 

Aero  and  Materials  Division 

Directorate  of  Systems  Engineering 

Deputy  for  Reentry  Systems 
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SECTION  1 


THE  PROBLEM 

As  an  ablating  nosetip  progresses  through  its  trajectory,  its  surface 
shape  changes  continuously.  Among  the  myriad  shapes  it  may  assume,  there  are 
certain  classes  of  shapes  for  which  the  flow  field  is  grossly  unsteady,  -that 
is,  the  flow  field  pulsates  with  a high  frequency.  Such  shape  induced 
flow  field  instabilities  may  seriously  compromise  the  realization  of  flight 
objectives,  possibly  even  resulting  in  structural  failure  of  the  nosetip.  This 
report  describes  the  results  of  an  experimental  and  analytical  effort  directed 
toward  improving  prediction  procedures  for  identifying  instability  onset  and, 
ultimately,  the  magnitude  and  frequency  of  the  pressure  fluctuations. 

Prior  to  this  study,  there  were  experimental  data  on  such  flow  field 
instabilities  for  idealized  shapes,  primarily  flat  faced  cylinders  and  high 
angle  cones  with  slender  spikes.  The  data  were  adequate  to  define,  at  least 
roughly,  the  envelope  of  instability  onset  for  the  idealized  shapes.  There 
was  no  such  data  base  fcr  shapes  characteristic  of  ablating  nosetips.  The 
objective  of  the  current  experimental  program  is  tc  fill  this  void  by  providing 
a data  base  on  which  a semi-empirical  model  could  be  constructed  defining  the 
envelope  of  instability  onset  for  shapes  of  interest  in  ablating  nosetip  tech- 
nology. An  additional  objective  is  to  provide  data  which  could  be  used  to 
correlate  pressure  pulsation  frequencies  and  magnitudes  between  ground  test  and 
flight  configurations. 

The  analytical  study  has  a two-fold  objective:  (1)  to  analytically 

identify  the  geometric  parameters  which  are  important  in  defining  instability 
onset  and  construct  a semi-empirical  correlation  for  the  envelope  of  instability 
onset,  and  (2)  to  determine  the  feasibility  of  constructing  a less  empirically 
dependent  flow  field  model  for  predicting  the  onset,  frequency,  and  magnitude 
of  pressure  fluctuations  resulting  from  gross  flow  field  pulsations. 

The  experimental  and  analytical  programs  are  discussed  in  Sections  2 and 
3,  respectively.  The  results  of  the  two  efforts  are  then  synthesized  into  a 
correlation  for  predicting  instability  onset  for  nosetip  configurations  (Section 
4).  The  report  is  concluded  with  a summary  of  the  important  results. 
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SECTION  2 


EXPERIMENTAL  PROGRAM 


2.1  BACKGROUND,  PROGRAM  OBJECTIVES,  AND  APPROACH 

As  already  noted,  almost  all  previous  experimental  efforts  in  this  area 

have  been  focused  on  idealized  shapes.  Particularly  noteworthy  are  the  studies 
1 2 

of  Maull  and  Wood  of  the  hypersonic  flow  over  spiked  cones  and  cylinders. 

Wood  defined  five  flow  regimes  depending  on  the  cone  half  angle  and  the  ratio 
formed  by  dividing  the  spike  length  by  the  model  base  diameter  (Figure  2-1) . 

Wood  hypothesized  that  for  very  blunt  spiked  bodies,  a necessary  con- 
dition for  instability  is  that  the  spike  protrude  beyond  the  position  that  the 
bow  shock  would  be  detached  if  there  were  no  spike.  This  hypothesis  is  consis- 
tent with  existing  experimental  data.  Maull  noted  the  dependence  of  the  enve- 
lope of  instability  onset  to  variations  in  Mach  number,  aspect  ratio,  and  model 
shoulder  radius.  He  also  noted  that  the  boundary  between  steady  and  oscillatory 
flow  correlates  with  the  sharp  cone  angle  for  shock  detachment.  On  the  basis 
of  a number  of  shadowgraphs  taken  during  the  tests,  Maull  constructed  a physi- 
cal model  of  the  pulsation  cycle  in  which  the  flow  field  instability  onset  is 
related  to  the  inability  of  a separated  shear  layer  to  mount  the  strong  adverse 
pressure  gradient  associated  with  the  strong  shock  created  at  the  shoulder  when 
the  conical  shock  detaches  from  the  conical  surface. 

The  geometries  tested  and  correlated  by  Maull,  Wood  and  most  others  were 
spiked  cylinders  and  large  angle  cones,  shapes  not  representative  of  configur- 
ations encountered  by  ablating  nosetips.  Prior  to  this  test  series  there  was 
a paucity  of  data  about  flow  instabilities  on  shapes  likely  to  occur  on  ablating 
nosetips . 

The  nature  of  these  tests  is  that  of  a screening  program  to  determine  the 
range  of  geometric  configurations  and  environmental  conditions  which  result  in 
large-scale  flow  field  pulsations  and  to  determine  the  frequency  and  magnitude 
of  the  resulting  pressure  oscillations.  This  required  modeling  the  ablated 
nosetip  shapes  so  that  they  are  characterized  by  am  acceptably  small  number  of 
geometric  parameters.  Then  low  temperature  ablator  data,  limited  shape  change 
predictions  for  flight  conditions,  and  Mach  number  effect  considerations  were 
used  to  determine  the  range  of  geometric  parameters  of  interest  in  this  test 
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program.  A total  of  38  different  shapes  was  tested.  Since  these  were  screen- 
ing tests , a few  shapes  were  selected  because  they  were  almost  certain  to  pul- 
sate or  to  remain  steady,  but  most  of  the  configurations  were  selected  because 
it  was  uncertain  whether  or  not  they  would  pulsate,  and,  therefore,  the  tests 
provide  the  needed  data  base  to  better  define  the  geometrical  configuration 
boundaries  corresponding  to  instability  onset. 

T'.»e  models  were  tested  in  the  Naval  Ordnance  Laboratory  tunnel  number  8 
in  Silver  Spring,  Maryland,  at  Mm  = 5.  The  tips  of  the  models  were  gritblast 
roughened  to  ensure  a turbulent  boundary  layer  on  the  trp  if  the  flow  remained 
attached.  Comparative  visual  and  tactile  inspection  of  the  roughened  tips 
indicate  that  the  roughness  was  nominally  in  the  x-2  mils  range.  Detailed 
roughness  characterization  studies  were  not  conducted. 

Each  model  was  instrumented  with  two  high  frequency  pressure  transducers. 
The  sting  was  instrumented  with  an  axial,  a yaw  plane,  and  a pitch  plane  accel- 
erometer. Pulsed  35  mm  still  Schlieren  (about  4/sec)  and  16  mm  high  speed 
shadowgraph  cinephotography  coverage  were  used  on  each  test,, 

2.2  MODEL  SELECTION 

It  is  not  easy  to  model  simply  the  surface  geometry  of  an  ablated  shape. 
However,  in  order  to  develop  and  use  a correlation  which  defines  the  geometri- 
cal dependence  of  instability  onset,  it  is  necessary  to  describe  the  surface 
geometry  with  relatively  few  parameters.  A careful  examination  of  low  tempera- 
cure  ablator  shapes  and  shape  change  prediction  results  led  to  the  conclusion 
that  the  surface  geometry  can  be  adequately  characterized  by  the  model  indicated 
in  Figure  2-2.  By  varying  the  parameters*  L/D,  Rn/D,  Rf/D,  Rs/D,  8^  and 
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Figure  2-2.  Schematic  of  Simplified  Model  of  Nosetip 


During  the  data  correlation  efforts  a different  definition  of  L was  determined 
to  better  correlate  the  data.  See  Section  2.6.3. 
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through  the  range  of  interest  determined  by  examining  the  low  temperature 

ablator  data  and  published  experiments  with  pulsating  flow  on  spiked  bodies , 

the  envelope  of  instability  onset  for  realistic  configurations  was  determined. 

With  such  a large  numbe,  'f  parameters,  it  was  unrealistic  to  test  all  possible 

combinations  of  interest.  Therefore,  nominal  values  of  R , L/D , R,/D,  R „/D,  0. , 

n t e i. 

and  02  were  determined.  The  tests  involved  varying  these  parameters  about  their 
nominal  values  to  determine  the  sensitivity  of  the  onset  and  intensity  of  flow 
field  instability  to  departures  from  the  nominal  values.  The  result  is  that 
the  test  objectives  were  achieved  by  testing  thirty-three  different  configur- 
ations with  the  following  range  of  parameter  variation: 
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L/D: 

0.1, 

0.2,  0.3,  0.4,  0 

Rn/D: 

* 

O 

• 

O 

0.15,  0.30 

Rf/D: 

0.0, 

0.05,  0.10,  0.20 

Rg/D: 

0.0, 

0.05,  0.10,  0.20 

The  test  matrix,  as  run,  is  given  in  Table  2-3,  Section  2.5. 

2.3  MODEL  DESIGN  AND  INSTRUMENTATION 

As  was  mentioned  in  Section  2.2,  the  nature  of  the  tests  described  here- 
in is  that  of  a screening  program  to  determine  the  range  of  geometric  configur- 
ations and  environmental  conditions  which  result  in  large  scale  flow  pulsations. 

A total  of  33  different  shapes  was  tested.  To  accomplish  the  test  objectives 
while  minimizing  model  expense  and  test  time  required  that  the  models  be  designed 
so  that  some  geometric  variations  could  be  realized  without  dismounting  the 
entire  model.  This  section,  2.3.1,  includes  a detailed  discussion  of  model 
design.  Sections  2.3.2  and  2.3.3  discuss  pressure  and  accelerometer  instrumen- 
tation and  finally,  photographic  coverage  and  data  acquisition  is  discussed  in 
Sections  2.3.4  and  2.3.5. 

2.3.1  Models 

In  order  to  minimize  model  expense  and  test  time,  the  models  were  construct- 
ed in  two  parts,  a tip  and  a transition  ring,  illustrated  in  Figure  2-3.  Table 
2-3  defines  the  pertinent  geometric  variables  for  the  various  transition  rings 
and  tips.  With  the  exception  of  transition  ring  number  8 and  tips  10  and  30, 
tips  are  characterized  by  four  geometric  variables:  R^,  R^,  L and  6^,  whereas 

transition  rings  are  characterized  by  Rs,  02  and  D (refer  to  Figure  2-2  or 
Table  2-3).  Tip  numbers  10  and  30  are  characterized  by:  R^j,  R^,  L,  0^  Rg,  and 
©2 . Transition  ring  number  8 was  merely  a cylindrical  sleeve.  Model  base 
diameter  (D)  was  6.0  inches  for  all  transition  rings.  The  tips  and  transition 
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rings  2 through  7 were  fabricated  from  a high  strength  carbon  steel  (Type  1215) ; 
transition  ring  8 was  fabricated  from  6061-T6  aluminum. 

Figure  2-4  illustrates  transition  ring  number  2,  and  Figures  2-5,  2-6, 
and  2-7  depict  tips  8,  6,  and  1,  respectively.  Some  of  the  important  design 
features  illustrated  in  these  figures  are: 

• Tolerance  between  tip  and  transition  ring  held  at  a 
nominal  .002  inch. 

• 0-ring  used  to  seal  interior  cavity  of  the  transition 
ring 

• Tips  held  in  place  via  4 set  screws  (refer  to  section  3-B, 

Figure  2-5) 

• Pressure  ports  for  6 static  pressure  taps  and  for  2 high  frequency 
pressure  transducers  were  machined  only  on  the  transition  ring. 

Obviously,  the  model  was  designed  to  allow  quick  changes  in  configuration. 

The  nosetip,  being  attached  via  four  set  screws,  could  be  changed  in  approximate- 
ly five  minutes.  The  transition  section,  housing,  two  high  frequency  pressure 
transducers  and  six  static  pressure  taps  could  be  changed  in  approximately 
one- half  hour. 


The  assembly  drawing  (Figure  2-8)  illustrates  the  model  wind  tunnel 
sting  interface,  which  includes,  from  left  to  right:  tip,  transition  ring, 

model  backplate,  model  holder,  adapter  plate,  instrumented  spacer,  and  sting. 
Static  pressure  tap  tubing  and  electrical  wires  from  the  Kulite  pressure  trans- 
ducers *.?ere  placed  along  the  center  of  parts  illustrated  in  Figure  2-8.  Elect- 
rical connectors  were  mounted  inside  the  transition  rings  so  that  wiring  could 
be  easily  connected  or  disconnected  to  the  high  frequency  pressure  transducers. 
The  capillary  tubes  from  the  static  pressure  taps  were  connected  inside  the 
model  holder  to  tubing  which  ran  through  the  sting  and  into  a bank  of  static 
pressure  transducers  located  outside  of  the  test  chamber. 


In  general,  the  static  and  high  frequency  pressure  transducer  ports  on 
the  face  of  the  transition  ring  arc  located  as  close  as  possible  to  the  flare 
so  that  if  there  is  substantial  boundary  layer  separation  the  measured  pressures 
will  be  in  the  separated  flow  region.  For  L = 2.4  inches  and  8j-  30  degrees, 
i.e.  tips  10  and  30,  it  was  not  practical  to  construct  the  models  in  two  parts 
while  still  maintaining  the  pressure  ports  on  the  transition  ring  and  locating 
the  pressure  ports  in  the  region  of  potential  separated  flow.  Therefore,  tips 
10  and  30  were  constructed  as  illustrated  in  Figure  2-9.  Figure  2-10  depicts 
Tip  27  Transition  Ring  8 interfaced  with  the  wind  tunnel  facility  sting. 
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To  insure  a turbulent  boundary  layer  by  the  first  shoulder  on  the  model 
for  a nominal  free  stream  unit  Reynolds  number  of  20  x 106/ft,  the  model  tips 
were  paLtially  roughened.  A nominal  peak  to  valley  roughness  height  of  1-2  mils 
was  chosen  based  on  the  transition  criteria  recommended  in  Reference  3.  The  1-2 
mils  surface  roughness  was  obtained  via  grit  blasting  the  model  tips. 

One  additional  point  needs  to  be  made  regarding  the  fit  between  the  tip 
and  transition  ring.  In  general,  the  surface  of  a model  at  the  point  where  the 
tip  and  transition  ring  meet  is  sligntly  discontinuous.  For  transition  rings 
2,5,  and  6,  a step  size  on  the  order  of  1 to  2 mils  resulted  from  a slight 
breaking  of  the  material  on  the  sharp  corner  on  the  transition  ring.  Note, 
however,  that  this  surface  discontinuity  is  of  the  same  order  or  less  than 
the  surface  roughness  and,  therefore,  should  have  had  no  adverse  effect  on  the 
results. 

2.3.2  Pressure  Instrumentation 

The  models  tested  during  this  program  were  instrumented  with  static 
pressure  taps  and  high  frequency  pressure  transducers. 

Static  Pressure  Instrumentation 


Each  transition  ring  (except  transition  ring  8)  was  fitted  with  6 pres- 
sure tubes:  two  on  the  face  served  as  reference  pressures  for  the  high  frequency 

transducers,  and  4 aft  of  the  shoulder.  In  addition,  one  static  pressure  tube 
was  placed  inside  the  model  to  monitor  the  internal  pressure. 

As  illustrated  in  Figure  2-4  (see  Section  2.3.1)  , the  static  pressure 
taps  on  the  face  of  the  transition  ring  were  constructed  by  drilling  a 0.055 
inch  diameter  hole  normal  to  the  surface  and  counterboring  a 0.12b  inch  diameter 
hole  to  within  1/4  of  an  inch  from  the  surface.  Stainless  steel  tubing  (outside 
diameter. 0. 125  inch  inside  diameter  0.095  inch)  was  then  brazed  into  the  0.128 
inch  diameter .hole.  The  static  pressure  tap  holes  on  the  cylindrical  section  of 
the  transition  ring  were  constructed  by  drilling  a 0.096  inch  diameter  hole  and 
brazing  stainless  steel  tubing  having  outside  and  inside  diameters  of  0.093  inch 
and  0.061  inch,  respectively.  In  the  model  holder  (c.f.  Figure  2-10)  the  stain- 
less steel  tubing  originating  on  the  cylindrical  section  of  the  transition  ring 
was  connected  to  0.125  inch  OD  stainless  steel  tubing. 

Table  2-1  defines  the  pressure  tap  locations  for  all  transition  rings. 


The  pressure  tap  holes  are  referred  to  as  ’1  through  P6.  PI  and  P2  are  located 
on  the  face  of  the  transition  ring,  in  the  yaw  plane,  whereas  P3  through  P6 
are  located  on  the  cylindrical  section  90  degrees  apart.  The  cylindrical 
coordinate  system  used  to  describe  the  pressure  tap  location  is  illustrated  in 
the  sketch  associated  with  Table  2-1.  For  all  transition  rings,  the  origin  of 
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• -neasured  frora  stagnation  point 

- strearvise  distance  to  these  pressure  ports  for  tips  10  and  30  is  2.797  inches  and  2.228  inches,  respectively 

• streanwise  distance  to  these  pressure  ports  for  tips  10  and  30  is  4.018  inches  and  3.757  inches,  respectively 
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the  coordinate  system  is  located  at  the  model  centerline  at  the  beginning  of 
the  transition  ring.  The  streamwise  or  axial  distance  from  the  model  stag- 
nation point  can  be  readily  computed  for  any  tip  transition  ring  configuration 
using  the  geometric  information  given  in  Table  2-3.  As  discussed  in  Section 
2.3.2,  model  tips  10  and  30  were  solid  models  and  incjude  both  static  and  high 
frequency  pressure  taps  on  the  fore-cone  of  tne  tip.  Table  2 -lb  defines  the 
location  of  these  pressure  taps  in  terms  of  radial  distance  from  the  model 
centerline  (r)  , meridian  angle  ( <J>)  , and  axial  distance  from  model  stagnation 
point.  Also  included  is  the  streamwise  distance  to  the  pressure  ports.  $ = o° 
corresponds  to  the  leeward  side  for  the  angle  of  attack  runs. 

Dynamic  Pressure  Instrumentation 

Two  high  frequency  pressure  transducers  were  located  on  the  transition 
ring  face  in  the  pitch  plane.  Table  2-1  gives  the  transducer  (labeled  K1  and 
K2)  locations  in  terms  of  the  coordinates  z,  r,  and  $•. 

Kulite  pressure  transducers , manuf acted  by  Kulite  Semiconductor  Products , 
Inc. , were  selected  for  obtaining  the  high  frequency  pressure  response  on  the 
models.  The  model  number  of  the  Kulite  high  frequency  pressure  transducers 
is  XTEL-190-100 . Figure  2-11  gives  the  nominal  specifications  of  the  trans- 
ducers . These  transducers  u-e  pie2oresistive;  four  piezoresistive  elements 
are  diffused  on  a silicone  diaphragm  forming  a four  arm  Wheatstone  bridge. 

The  transducers  have  a rated  pressure  of  100  psi  and  a maximum  pressure  of 
200  psi.  The  Kulite  transducers  arc  equipped  with  a one  inch  reference  tube 
and,  therefore,  are  classified  as  pressure  gauges.  As  mentioned  above,  a 
static  pressure  tube  was  positioned  in-side  each  model  so  that  the  interior 
pressure  could  be  measured,  thus  allowing  one  to  obtain  a relative  comparison 
between  pressures  measured  via  Kulite  gauges  and  capillary  tubes. 

The  operating  temperature  for  the  XTEL-90-100  transducer  is  0°F  to  500°F. 
These  transducers  are  equipped  with  a thermal  compensation  unit  providing 
temperature  compensation  from  80°F  to  450°F.  The  electric  circuit  and  external 
dimensions  of  the  XTEL-90-100,  including  the  thermal  compensation  unit,  is 
illustrated  in  Figure  2-11. 

The  Kulite  pressure  transducers  were  mounted  in  a ported  configuration 
(diameter  of  port  = 0.060  inch,  length  of  port  = 0.10  inch),  as  illustrated  in 
Figures  2-4  and  2-11.  Kulite  transducers  were  dynamically  tested  in  a shock 
tube  to  determine  the  transducer  response  to  a step  input  in  pressure  (approxi- 
mately 30  psi).  The  results  of  these  tests  indicate  that  a Kulite  gauge,  mounted 
in  the  ported  configuration  given  above,  experienced  output  oscillations  with  a 
frequency  of  approximately  13  KHz,  with  the  oscillations  being  naturally  damped 
in  300  microseconds.  Similarly,  a Kulite  gauge  mounted  flush  in  the  shock  tube 
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BRIDGE  TYPE 

RATED  PRESSURE 
MAXIMUM  PRESSURE 
SENSITIVITY-NOMINAL 
EXCITATION 
INPUT  IMPEDANCE 
OUTPUT  IMPEDANCE 
ZERO  BALANCE 

COMBINED  NON-LINEARITY  AND 
HYSTERESIS 

REPEATABILITY 

COMPENSATED  TEMPERATURE  RANGE 

OPERATING  TEMPERATURE  RANGE 

CHANGE  OF  SENSITIVITY  WITH 
TEMPERATURE 

.CHANGE  OF  NO-LOAD  OUT  ’JT  WITH 
TEMPERATURE 

NATURAL  FREQUENCY  . 

ACCELERATION  SENSITIVITY 


Fully  active  four  am  Wheatstone  bridge 
diffused  on  silicon  dlaphragn 

100  psl 

200  psl * 

15  nv/psl 

1U  vDC  (20vDC  Maxlmtan) 

2500  ohms  (non.) 

MOO  ohms  (non.) 

± 5X  (rated  pressure) 

IX  (rated  pressure) 

0.5X  ^ 

80°F  to  450°F  ' 

0°F  to  500*F 

± 5X  within  the  compensated  temperature  range 

± 10X  of  rated  pressure  within  the  compensated 
temperature  range 

Greater  than  100  KHz 

Perpendicular:  0.0005  psi/g 
Transverse:  0.00005  psi/g. 


‘Thermal  set  is  defined  as  the  change  in  specification  doe  to  • the;  msl  cycle  from  the  km  end  to  the  high 
end  of  the  corr.penuted  temperature  range. 
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results  in  output  oscillations  on  the  order  of  the  natural  frequency  of  the 
gauge,  100  KHz.  Output  oscillations  damped  out  after  600  microseconds  in  the 
flush  mounted  configuration. 


High  Frequency  Pressure  Instrumentation 


The  Kulite  pressure  transducers  used  during  these  tests  were  statically 
calibrated  by  the  manufacturer  and  Aerotherm.  Both  pre  and  post  test  calibrations 
were  statically  performed  at  the  tunnel  facility  upon  installation  in  a model. 


Table  2-2  summarizes  the  results  of  the  offsite  calibrations  for  transducers 


D46,  D47,  D48 , D19,  and  D50.  The  transducer  sensitivity  is  given  in  terms  of 
Mv/psi  per  bridge  excitation  voltage.  These  transducers  were  calibrated  with 
a bridge  excitation  voltage  of  DOvDC.  Also  included  in  Table  2-2  sure  (1)  tlje 
specific  transducers  which  were  mounted  in  the  7 transition  rings  and  two  solid 
nosetips,  (2)  transducer  labels  associated  with  mounting  (K1  and  K2  refer  to  0 
and  180  degree  meridians,  respectively),  (3)  zero  offset  of  the  transducers, 

(4)  estimation  of  the  thermal  effect  on  transducer  sensitivity  and  (5)  thermal 
effect  on  zero  pressure  output  variation. 


In  addition  to  the  two  Kulite  transducers  mounted  on  the  model,  a Kulite 
gauge  (K3)  was  flush-mounted  on  the  nozzle  wall  to  obtain  pressure  fluctuations 
of  the  nozzle  side  wall  boundary  layer.  The  data  obtained  from  K3  was  to  be 
reduced  by  NOL  personnel  in  order  to  quantify  the  magnitude  of  noise  generated 
by  the  turbulent  side  wall  boundary  layer. 


All  three  Kulite  transducers  were  attached  to  a facility  D.C.  power 
supply;  the  nominal  input  voltage  to  the  transducer  was  6.95  volts.  Preston 
model  8300  XWBRC  amplifiers  were  used  with  the  high  frequency  pressure  trans- 
ducers. The  band  width  on  the  amplifiers  was  set  at  the  maximum,  100  KHz,  and 
the  gain  was  adjusted  to  give  an  appropriate  input  signal  to  the  FM  tape 
recorder  (the  tape  recorder  is  described  in  Section  2.3.5). 


2.3.3  Accelerometer  Instrumentation 


Three  Endevco  2220C  high  frequency  accelerometers  were  mounted  aft  of  a 
model  in  the  instrumented  spacer,  as  shown  in  Figures  2-8  and  2-10.  A schematic 
of  the  instrumented  spacer  is  given  in  Figure  2-12.  Accelerometers  1,  2 and 
3 measured  accelerations  in  the  vertical,  horizontal  and  axial  directions, 
respectively.  The  Endevco  2220C  accelerometers  have  a range  from  0 to  5000g 
and  a natural  frequency  of  50  KHz. 


Endevco  2710  charge  amplifiers  were  used  to  amplify  the  output  of  the 
piezoelectric  accelerometers.  This  Endevco  charge  amplifier  has  a frequency 
range  of  2 to  20,000  Hz,  with  a maximum  gain  of  10  volts  per  pico  coulomb 
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j input.  For  all  uses  in  this  test  series,  the  charge  amplifiers  were  set  at  the 

] maximum  frequency  response  and  approximate xy  one  half  of  maximum  gain. 
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2.3.4  Photographic  Coverage 

Both  35  mm  pulsed  still  Schlieren  and  16  mm  high  speed  shadowgraph  movie 

1 

coverage  were  used  in  order  *-.o  visualize  large  scale  flow  oscillations  and  to 
| determine  whether  boundary  layer  separation  existed  when  the  flow  was  steady. 

Figure  2-13  schematically  illr  trates  the  camera  setup.  As  shown,  two  35  mm 
pulse  cameras  were  used,  one  viewing  the  model  pitch  plane  through  the  off-axis 
•*Z*  Schlieren  system,*  the  other  photographing  a millisecond  counter  located 
outside  of  the  tunnel  test  section.  Both  35  mm  cameras  are  Flight  Research 
(Model  IV)  multi-data  cameras  which  were  equipped  with  a 9 inch  len3.  The  35mm 
camera  viewing  the  model  through  the  Schlieren  system  was  synchronized  to  a 
pulsed  strobe  light  (General  Radio,  model  1535)  which  was  interfaced  with  the 
Schlieren  system.  The  35  mm  timing  camera  was  also  synchronized  to  the  35  mm 
viewing  camera.  Both  35  mm  cameras  were  run  at  approximately  4 frames  per 
second. 

i 

! 

The  high  speed  camera  was  a 16  mm  Hycam  movie  camera  which  was  run  at 
20,000  frames  per  second.  The  Hycam  camera  was  operated  in  a half  frame  mode 
resulting  in  a frame  size  of  0.15  inch  by  0.41  inch.  This  camera  was  equipped 
with  a j inch  lens  and  was  focused  through  the  Schlieren  system  and  had  the 
same  view  of  the  model  as  the  35  mm  pulse  camera.  The  lighting  for  the  Hycam 
camera  was  obtained  from  the  Schlieren  system,  i.e.  from  the  high  pressure 
mercury  vapor  lamp,  but  the  lens  system  was  oriented  to  yield  shadowgraph  type 
coverage. 

I In  general,  35  mm  pictures  were  taken  during  an  entire  model  test,  which 

was  typically  about  10  seconds . A problem  regarding  proper  synchronizing  of 

j the  35  mm  pulse  cameras  and  the  strobe  light  was  e untered  early  in  the  test 

matrix.  This  problem  resulted  in  varying  frame  exposures  and  incorrect  filming 
of  the  millisecond  counter.  Despite  varying  frame  exposures,  the  35  mm  photo- 
graphs of  the  model  and  flow  field  were  of  sufficient  quality  to  accomplish 
the  test  objectives.  Moreover,  an  accurate  temporal  history  of  the  35  mm 
photographs  was  not  a necessity  in  terms  of  data  reduction,  and,  therefore, 
the  incorrect  filming  of  the  millisecond  counter  did  not  compromise  realization 
of  the  program  objectives. 

High  speed  movie  films  were  obtained  on  all  runs,  except  angle  of  attack 
runs  and  variable  Reynolds  number  runs  for  which  the  freestream  unit  Reynolds 
number  was  less  than  20  x 10*/ft.  The  Hycam  uamera  was  started  shortly  after 
model  insertion  in  the  test  stream,  and  coverage  lasted  approximately  1.5  seconds 
thereafter. 
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2.3.5  Data  Acquisition 


As  can  be  inferred  from  the  above  sections,  recorded  data  for  the  Series 
G wind  tunnel  tests  included: 

• Static  pressure  measurements 

• High  frequency  pressure  measurements 

• Accelerometer  measurements 

• 35  mm  film  coverage 

• High  speed  (20,000  frames/sec)  16mm  film  coverage 

A brief  description  of  the  data  recording  system  is  presented  in  the  following 
paragraphs . 

At  the  NOL  wind  tunnel  facility,  a 14  channel  digital  recorder  is  used 
for  recording  the  test  time,  supply  pressure,  supply  temperature  and  static 
pressures.  This  recorder  has  a sampling  or  sweep  rate  of  0.056  second.  Test 
time,  supply  temperature,  supply  pressure,  and  model  position  are  recorded  on 
separate  channels,  leaving  10  channels  available  for  recording  static  pressures. 

Static  pressure  tubing  from  the  model  was  attached  to  NOL  static  pressure 
transducers,  which  were  monitored  by  the  digital  recorder.  Pressure  taps  on 
the  aft  cone  of  the  model  were  attached  to  transducers  with  a pressure  range 
from  0 to  50  p3i,  whereas  pressure  taps  on  the  cylindrical  portion  of  the  model 
were  attached  to  transducers  with  a pressure  range  from  0 to  5 psi.  The  static 
pressure  transducers  were  calibrated  every  2 to  3 runs. 

The  digital  recording  system  places  the  data  on  magnetic  tape,  which  was 
subsequently  reduced  according  to  the  calibrations.  During  each  test,  digital 
displacement  of  pressure  transducer  emf  output  as  a function  of  time  was  available 
for  four  channels. 

High  frequency  Kulite  pressure  transducer  and  Endevco  accelerometer  out- 
put was  amplified  and  recorded  on  a 14  channel  FR-1300  Ampex  tape  recorder. 

This  tape  recorder  was  run  in  a FM  mode  with  a tape  drive  speed  of  60  inches/sec- 
ond  which  results  in  a flat  frequency  response  from  0 to  20,000  cycles/second. 

In  addition  to  six  channels  of  tl  e tape  recorder  taken  by  three  Kulite 
pressure  transducers  and  three  Endevco  accelerometers,  an  IRIG  B signal  and  a 
50  kHz  sine  wave  were  recorded  on  two  of  the  remaining  8 channels . The  IRIG  B 
timer  signal  generator  was  manufactured  by  Datatron,  Inc.,  Model  3000.  The 
sine  wave  was  produced  by  a Hewlett-Packard  wave  form  generator.  An  oscilli- 
scope  and  oscillograph  were  also  available  to  visualize  or  play  back  pressure 
transducer  or  accelerometer  response. 
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2.4  WIND  TUNNEL 

The  Series  G wind  tunnel  tests  were  conducted  in  NOL  Tunnel  Number  8, 
which  io  an  intermittent  blowdown  hypersonic  wind  tunnel.  For  this  tunnel, 
air  is  the  working  gas  and  is  compressed  and  stored  in  high  pressure  containers. 
In  transit  to  the  nozzle,  the  air  is  heated  to  the  desired  supply  temperature 
by  passing  it  through  a pebble  bed  heater.  The  tunnel  has  an  open  jet  test 
section  and  a constant  area  diffuser.  The  Mach  number  for  this  test  series 
was  Mm  = 5.  The  associated  nozzle  is  of  rectangular  cross  section  with  throat 
dimensions  of  0.6398  inch  by  16.00  inches,  ana  exit  dimensions  of  17.37  inches 
by  16.97  inches. 

At  a Mach  number  of  5,  the  Reynolds  number  operating  envelope  for  Tunnel 
Number  8 is  presented  in  Figure  2-14  as  a function  of  supply  temperature  (T  ) 
and  supply  pressure  (Po) . Superimposed  on  the  figure  are  the  four  nominal  oper- 
ation conditions  which  were  run  during  the  test  series.  The  lowest  Reynolds 
number  runs  were  not  achieved  because  of  breakdown  of  the  tunnel  flow.  T“e 
actual  supply  conditions  are  given  in  Table  2-3  (see  Section  2.5^  . 

Models  were  aligned  in  the  tunnel  test  stream  by  manually  checking  the 
anqularity  of  the  models  in  the  pitch  and  yaw  plane  with  a clinometer. 

2 . 5 TEST  MATRIX 

The  preliminary  data  analyses  indicated  that  6^,  02,  and  L/D  are  the 
most  important  parameters  (verified  in  these  tests) , and  the  test  matrix  empha- 
sizes variations  in  those  three  variables.  Effects  cf  Rn/D»  Rf/D,  and  Rg/D  were 
then  determined  for  nominal  values  of  0^,  02,  and  L/D. 

One  may  argue  that  the  important  nose  radius  parameter  is  R rather  than 
Rn/D  since  it  is  likely  that  the  value  of  Rn  corresponding  to  a stable  shape  is 
independent  of  the  represen tative  nosetip  diameter  D.  Consideration  of  boundary 
layer  transition  might  lead  one  to  the  same  conclusion  (Section  2.6.3).  How- 
ever, geometric  similarity  considerations  might  lead  one  to  conclude  that  the 
ratio  Rn/D  is  more  important  with  respect  to  instability  onset.  Nose  radius  is 
an  important  factor,  but  the  precise  reason (s)  is  not  well  understood. 

The  tests  were  run  at  nominal  total  temperatures  of  about  300°F  and 
a free  stream  unit  Reynolds  number  of  20  x 10*  per  foot.  The  effect  of  vary- 
ing Reynolds  number  was  studied  on  three  geometries  by  decreasing  the  free 
stream  unit  Reynolds  number  to  10  and  5 x 10*  per  foot.  It  was  initially  intend- 
ed that  data  also  be  obtained  at  Re^/foot  = 2.5  x 10s,  but  the  tunnel  flow  broke 
at  that  low  total  pressure.  Most  configurations  were  tested  only  at  zero  degrees 
incidence,  but  three  were  also  tested  at  a = 2 1/2°  and  5C  incidence  to  provide 
data  on  the  effect  of  angle  of  attack  on  flow  instability. 
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The  test  matrix,  as  run,  is  given  in  Table  2-3.  Also  included  in  Table 
2-3  is  a column  summarizing  the  type  of  flow  field/pressure  response  observed 
during  the  test.  Those  results  are  discussed  in  Section  2.6.2. 

2 . 6 RESULTS 

The  results  of  the  Series  G tests  ?-e  discussed  in  this  section.  The  dis- 
cussion begins  with  an  overview  of  the  test  results  in  terms  of  the  basic  find- 
ings of  the  experiments.  This  is  followed  by  a discussion  of  the  different 
types  of  flowfields  observed  in  the  motion  pictures  and  still  photos  and  how 
these  relate  to  the  mean  and  fluctuating  pressure  data.  Preliminary  attempts 
to  correlate  the  data  are  summarized.  Further  discussions  >f  the  correlation 
efforts  and  how  they  are  directly  affected  by  theoretical  considerations  are 
given  in  Sections  3 and  4. 

2.6.1  Overview  of  Findings 

Important  objectives  of  the  experiments  are  to: 

• Identify  critical  parameters  which  influence  flow  stability 

• Determine  effects  of  Mach  number,  Reynolds  number,  and 
geometry  upon  stability  limits 

• Identify  ways  in  which  theory  can  be  used  to  extrapolate  to 
new  conditions 

• Determine  best  means  of  correlating  present  results  and  other 
existing  data  for  application  to  ablating  nosetip  technology. 


Table  2-3  presents  a summary  of  the  test  configurations,  test  conditions, 
and  results  for  these  experiments . Included  in  this  table  is  a note  relating 
to  the  type  of  flow  observed  at  the  test  conditions  of  interest.  Four  such 
types  of  flow  are  noted:  completely  steady  (CS) , steady  separated  (SS) , fully 

pulsating  (FP) , and  metastable  (MS) . * jomplete  description  of  each  type  of 

flow  is  discussed  in  Section  2.6.2;  t}  i pulsed  35  mm  Schlie  an  photographs 
are  shown  in  Figure  2-13.  These  flow  ssifications  were  determined  by  ob- 
serving, in  detail,  the  16  mm  high  speed  films,  the  35  mm  photos  and  the  high 
frequency  pressure  and  accelerometer  data.  The  large  body  of  information  con- 
tained in  the  16  mm  films  was  reduced  by  recording  a tunning  verbal  commentary 


on  a tape  recorder  and  simultaneously  sketching  the  significant  observations. 
The  recording  was  later  transcribed  into  a series  of  notes  which  could  be 
easily  referred  to  in  order  to  recall  the  character  of  the  flow. 


Completed  Steady  (L/D  = .2)  Fully  Pulsating  (L/D 


The  four  types  of  flows  are  visually  very  distinct  from  each  other 
(Figure  2-15) . The  frames  shown  illustrate  many  of  the  characteristic  details 
of  the  flows  observed.  In  addition,  the  character  of  the  pressure  fluctuations 
felt  on  the  flare  face  are  indicated  by  the  insert  in  the  lower  left  hand 
corner  of  the  frames  in  Figure  2-15. 

Before  examining  the  detailed  results  presented  ir.  the  next  section, 
the  following  points  represent  a summary  of  the  general  findings  of  this  study. 

• The  shape  variables  6^,  02,  and  1/D  are  of  primary  importance  in 
determining  the  limits  of  stability. 

• The  variable  0^  has  special  significance  in  determining  the 
stability  limits. 

• The  fundamental  frequency  of  the  strong  oscillations  present  was 
about  1100  cps  which  corresponds  to  a Strouhal  number  of  about 
.20  which  agrees  with  other  data  in  the  literature. 

• Many  of  the  M3  type  flows  represent  the  most  severe  cases  of  pressure 
oscillations  on  the  nosetip  with  the  maximum  pressure  exceeding  the 
mean  by  more  than  a factor  of  two. 

• On  an  empirical/theoretical  basis  it  was  concluded  that  the 
stability  limits  are  not  strongly  affected  by  free  stream  Mach 
number  above  M„  = 5. 

CO 

• The  effects  of  Reynolds  number  and  tip  radius  have  some  influence  on 
stability,  but  probably  of  a secondary  nature. 

• The  effects  of  flare  radius,  R^,  and  shoulder  radius,  R£,  upon 
stability  limits  are  weak  or  negligible  for  the  range  of  these 
variables  tested. 

• The  effects  of  angle  of  attack  are  ak  or  negligible  for  the  range 
of  parameters  tested. 


2.6.2  Classification  of  Flows  Observed 

As  mentioned  previously  four  basic  types  of  flow  were  observed  in  these 
experiments : 

• Steady  separated,  SS 

• Completely  steady,  CS 

• Fully  pulsating,  FP 

• Metastaile,  MS 
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These  four  basic  types  of  flow  possess  radically  different  behavior  in 
terms  of 

• Visual  steadiness 

• Magnitude  of  pressure  oscillations 

• Magnitude  of  accelerometer  excitation 

Each  of  these  four  types  of  flow  are  discussed  in  the  following  paragraphs. 

2. 6. 2.1  Completely  Steady  (CS) 

Completely  steady  flow  means  that  the  shock  waves,  high  frequency 
pressure  transducer  response  and  accelerometer  response  are  all  completely 
steady.  Four  different  runs  exhibiting  completely  steady  flow  behavior  are 
shown  in  Figure  2-16.  These  flows  can  be  subclassed  depending  on  whether  the 
model  tip  protrudes  or  does  not  protrude  through  the  main  body  bow  shock. 

All  cases  in  which  the  tip  does  not  protiude  through  the  bow  shock  are  expected 
to  be  completely  steady.  A typical  case  is  run  41,  Figure  2-16a. 

Usually  when  the  model  tip  protudes  through  the  main  body  bow  shock, 
either  the  boundary  layer  separates  or  the  flow  pulsates.  However,  there 
were  three  tests  in  which  the  boundary  layer  remained  attached  and  the  flow 
field  was  apparently  completely  steady*.  Examples  are  runs  13,  39,  and  40, 
Figure  2-16. 

2.6. 2.2  Steady  Separated  (SS) 

The  photo  shown  in  Figure  2-17  taken  from  Run  20  is  typical  of  a steady 
separated  flow.  The  separation  on  the  tip  for  this  type  invariably  occurs  at 
a point  approximately  tangent  with  a line  drawn  from  the  tip  to  the  shoulder 
! of  the  model.  The  photograph  also  shows  a demarcation  line  between  the  flow 

I in  the  shear  layer  and  the  flow  in  the  separation  bubble  owing  to  the  density 

i 

gradient  that  exists  between  these  two  zones.  As  will  be  pointed  out  in  a 
later  section  this  demarcation  line  position  agrees  vrith  a prediction  of  the 
dividing  streamline  based  upon  the  shock  inclination  angle  off  the  tip.  For 
all  cases  of  steady  separated  flows  observed  the  reattachment  occurs  at  or 
very  near  th»  shoulder,  which  agrees  with  the  separation  model  discussed  in 
Section  3.0.  Occasionally,  for  many  of  the  steady  separated  flows  there 
appeared  a very  short  burst  of  unsteady  flow;  however,  this  only  lasted  for 


•Instrumentation  was  not  adequate  to  determine  very  small  localized  separation 
at  the  flare. 
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2-5  frames  of  the  high  speed  movies  and  the  flow  quickly  recovered,  regaining 
its  steady  appearance. 

2. 6. 2. 3 Fully  Pulsating.  FP 

The  basic  character  of  fully  pulsating  flow  is  depicted  in  Figure  2-18 
which  shows  a sequence  of  35  mm  photos  takers  during  Run  2.  A typical  sequence 
of  photos  was  chosen  which  corresponds  essentially  to  successive  frames  observ- 
ed in  the  16  mm  nigh  speed  films.  Fully  pulsating  flow  is  typified  by  continu- 
ous oscillations  and  no  indication  of  the  flow  "trying  to  become  steady".  The 
first  picture,  Figure  2- 18a,  shows  typically  the  beginning  of  the  oscillation 
cycle.  The  first  shock  >ff  the  tip  has  separated  flow  behind  it  as  measured 
by  the  local  instantaneous  shock  angle.  Such  estimates  are  only  approximate 
since  the  oblique  shock  relations  apply  truly  in  a moving  shock  frame  of 
reference.  The  second  frame  shows  the  shock  off  the  tip  increasing  in  inclina- 
tion due  to  an  expanded  growth  of  the  separated  region.  The  third  frame. 

Figure  2-18c,  shows  the  shock  off  the  tip  becoming  a strong  shock  as  the  separ- 
ation bubble  continues  to  grow.  Meanwhile,  the  secondary  shock  off  the  flare 
has  also  moved  forward  as  a result  of  the  bubble  growth.  Eventually,  (Figure 
2-13d)  the  first  and  second  shocks  nearly  merged  as  they  billow  off  in  front  of 
the  body.  At  this  point  the  accumulated  mass  in  the  bubble  begins  to  escape 
by  flowing  around  the  shoulder  of  the  body.  As  this  relief  occurs,  due  to  the 
shrinking  of  the  bubble  the  shock  structure  begins  to  collapse.  Figure  2-18e. 
Eventually,  the  shock  structure  returns  to  that  shown  in  Figure  2-18a  and  the 
sequence  of  events  is  repeated. 

2.6. 2. 4 Metastable  (MS) 

The  sequence  of  three  photos  shown  in  Figure  2-19  illustrates  a metasta- 
ble or  separated/fully  pulsating  flow  taken  from  Run  3°.  This  particular  type 
of  flow  is  characterized  by  having  alternating  periods  of  steady  separated  flow 
and  fully  pulsating  flow.  The  periods  of  steady  separated  flow  are  generally 
longer  in  duration.  This  type  of  behavior  is  most  difficult  to  explain  on  a 
theoretical  or  phenomenological  basis.  Figure  2-19a  shows  the  flow  pattern 
during  a period  of  nearly  steady  separated  flow.  At  some  later  instant  of 
time,  a ripple  in  the  shock  off  the  tip  appears,  indicating  that  some  sort  of 
disturbance  is  present.  This  disturbance  then  grows  in  intensity  throughout 
the  shear  layer  running  from  the  tip  to  the  shoulder.  Eventually  all  semblance 
of  a steady  separation  breaks  down  as  shown  in  Figure  2-19b.  This  flow  pattern 
resembles  the  start  of  a full  pulsation  sequence  and,  indeed,  the  oscillatory 
flow  pattern  which  ensues  is  fully  pulsating  in  nature  as  shown  by  another 
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typical  photo  from  the  oscillatory  mode  in  Figure  2-19c.  These  metastable 
flows  are  particularly  useful  in  defining  the  boundaries  between  fully  pulsating 
and  steady  flow. 

The  preceeding  remarks  describing  the  four  basic  types  of  flow  deal 
primarily  with  the  interpretation  of  observations  that  were  made  with  the  high 
speed  movies  and  the  still  photos.  A later  section  will  discuss  some  of  the 
quantitative  data  obtained  and  how  it  relates  to  the  visually  observed  flows. 


2.6.3  Effect  of  Geometric  and  Flow  Variables  Upon  Stability 


As  mentioned  previously,  a major  objective  of  the  present  test  series 
is  to  investigate  the  effects  of  several  geometric  and  flow  variables  upon  the 
stability  of  the  flow.  The  effect  of  the  various  parameters  investigated  is 
summarized  below. 

• 1/D  - the  results  of  these  tests  indicated  that  for  small  values  of 
1/D,  especially  if  the  tip  protrudes  not  at  all  or  only  slightly 
through  the  main  blunt  body  shock,  the  flow  is  likely  to  be  compiet'.ly 
steady  (CS) . On  the  other  hand  if  the  value  of  1/D  is  sufficiently 
large,  then  the  flow  is  likely  to  be  steady  separated  (SS) . However, 
if  the  proper  combination  of  0^  and  02  exist  and  if  1/D  falls  within 
the  limits  described  earlier  the  flow  is  likely  to  be  metestable  (MS) 
or  fully  pulsating  (FP) . 

• ®2  “ based  upon  theoretical  considerations  (see  Section  3.0)  if 
the  flare  angle  is  less  than  the  critical  angle  for  shock  detach- 
ment, the  flow  is  stable.  If  the  flare  angle  is  greater  than  the 
critical  angle  the  flow  may  be  unsteady  depending  upon  the  combin- 
ation of  0^  and  1/D  present. 

• 0 - apparently  if  0,  is  sufficiently  large,  the  flow  is  stable 

* - 1 I 

regardless  of  1/D  and  62-  Obviously,  if  = cot  (2^)  the  body 
is  a single  cone/cylinder  and  the  flow  is  steady.  Conversely,  the 
flow  has  increasing  tendency  to  be  unsteady  if  0^  is  small  and  1/D 
and  02  fall  within  the  limits  mentioned. 

• Rn/°  ~ there  is  a tendency  for  the  stability  of  the  flow  to  increase 
with  increasing  values  of  R^/D.  Figures  2-17  and  2-19  illustrate 
that  for  Runs  19  and  20  the  ability  of  the  flow  to  remain  steady  de- 
finitely increased  when  the  value  of  R _/D  was  increased  from  0.0  to 

05.  A possible  explanation  might  be  that  for  small  values  of  Rn/D 
the  boundary  layer  and,  hence,  the  shear  layer  is  laminar  just  after 
separation,  whereas  for  larger  values  of  Rn/D  the  boundary  layer 
and  shear  layer  are  turbulent  just  after  separation. 
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In  the  former  case,  transition  would  occur  along  the  shear  layers, 
resulting  in  a greater  tendency  toward  unsteadiness^. 

• R^/D  - the  effect  of  the  variable  R^/D  was  either  weak  or  negligible 
for  the  range  of  this  variable  tested.  There  was  no  noticeable 
effect  as  R^/D  was  varied  from  0.0  to  0.20. 
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R /D  - the  effect  of  the  variable  R /D  was  also  either  weak  or 

S 8 

negligible  for  the  range  of  this  variable  tested.  There  was  no 

noticeable  effect  as  R /D  was  varied  from  0.0  to  0.20. 

s 

R /ft  - the  free  stream  R /ft  was  varied  from  5 to  20  x 10s.  At 
e e'  % 

Rfi/ft  * 2.5  x 10*  and  at  times  for  5.0  x 10  a steady  tunnel  flow 

could  not  be  established.  Of  particular  interest  is  the  series  of 

runs  from  Run  27  through  Run  31.  Here  the  Re/ft  was  continuously 

reduced  holding  other  variables  constant,  resulting  in  an  increasing 

tendency  for  the  flow  to  become  unsteady.  Figure  2-20.  Figure  2-20a 

shows  a frame  taken  from  Run  27  with  R /ft  - 20  x 10s.  The  flow 

e 

is  essentially  steady  separated.  Figures  2-20b  and  2-20c  show  two 
frames  taken  from  Run  28  for  which  Rg/ft  = 10  x 10s;  here  the  flow 
is  shown  to  be  metastable.  Figures  2-20d  and  2-20e  show  two  frames 
taken  from  Run  30  where  Rg/ft  « 5 x 10s,  resulting  in  a fully  pulsat- 
ing flow.  Comparing  two  type  MS  flows  in  Table  2-3,  it  can  be  seen 
that  for  Run  28  the  flow  was  type  MS  at  a value  of  1/D  = .526  and 
Re/ft  = 10  x 10s  while  for  Run  32  the  flow  was  type  MS  at  a value 
of  1/D  = .367  at  R /ft  = 20  x 10s.  Superficially,  this  might  seem 
to  imply  that  for  geometry  considered, 

“^stable  , ilVitllj-“’ 

However,  further  inspection  of  the  geometry  of  Run  32  reveals  that 
oef£  = cot"1  (2  jjj) , where  ae££  is  the  effective  conical  body  created 
by  a separation  from  the  tip  to  the  shoulder  and  is  very  nearly  equal 
to  the  critical  angle  for  shock  detachment.  It  was  observed  that  MS 
type  flow  was  nearly  always  observed  under  such  a condition.  There- 
fore, for  Run  32  the  effect  of  Reynolds  number  is  obscurred  by  this 
geometric  condition,  and,  therefore,  the  coefficient  shown  in  equation 
2-1  is  probably  not  valid.  Therefore,  the  conclusion  is  qualitative 
and  indicates  that  flow  stability  increases  with  increasing  Reynolds 
number. 
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• a - a limited  number  of  angle  of  attack  variations  were  obtained 
during  Runs  25,  26,  42,  and  43 s here  the  angle  of  attack  was 
variously  set  at  0°,  2.5*,  and  5°.  Due  to  the  limited  amount  of 
angle  of  attack  testing,  the  results  are  somewhat  inconclusive. 
However,  Runs  25  and  26  did  seem  to  indicate  a trend  towards 
greater  instability  as  a increases. 

2.6.4  Anomolous  Types  of  Flow  Behavior 
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In  Idle  course  of  observing  the  high  speed  movies  and  still  photos  taken 
during  this  test  series  seve'  > instances  of  anomolous  flow  behavior  were 
noted,  falling  somewhat  out  c she  catagories  already  mentioned. 

The  metastable  type  flow  has  already  been  discussed  in  Section  2. 6. 2. 4. 
Figure  2-21  shows  other  examples  of  this  type  of  flow  behavior  taken  from  Runs 
6 and  38.  Owing  to  the  geometry  involved,  the  steady  shock  pattern  is  consider- 
ably different  between  the  two  runs,  and  Run  38  appears  to  have  a "blended* 
appearance  rather  than  the  two  shock  structure  of  Run  6.  The  other  point  of 
interest  is  that  during  the  oscillation  the  movement  of  the  strong  shock  for 
Run  38  is  considerably  less  chan  that  for  Run  6.  As  will  be  discussed  later, 
this  is  apparently  correlative  to  the  magnitude  of  the  variation  of  pressure 
on  the  flare  face. 

The  photos  shown  in  Figure  2-22  were  taken  from  Run  24  which  exhibited  a 
unique  flow  behavior.  These  results  indicate  that  although  the  flow  was 
essentially  of  the  FP  type  there  were  short  periods  of  time  when  the  flow  stead- 
ied out  and  appeared  as  shown  in  Figure  2-22a.  That  is,  these  results  are 
somewhat  similar  to  Run  38  except  that  the  periods  of  steady  flow  were  short. 
Note  also  that  when  the  flow  was  steady  the  tip  shock  was  not  "blended"  as  in 
Run  38  but  instead,  appears  as  two  distinct  shocks.  Thus,  Run  24  might  be  also 
classified  as  metastable,  or  perhaps  even  meta-unstable! 

Two  other  flow  anomolies  worth  mentioning  are  shown  in  Figure  2-23. 

Figure  2-23a  shows  a photo  of  Run  7 which  was  a metastable  type  flow.  Here  a 
disturbance  in  the  shear  layer  is  clearly  shown  to  begin  growing.  This  eventu- 
ally results  in  a period  of  oscillations.  Figure  2-23b  shows  a very  interesting 
shock  structure  which  occurred  during  a fully  pulsating  flow.  Run  9.  Note  the 
abrupt  cone-like  bulge  in  the  shock  structure  extending  ahead  of  the  main  oscil- 
lating shock  as  shown  in  the  lower  half  of  the  photo.  No  explanation  has  been 
* found  for  this  unusual  behavior. 
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2.6.5  Quantitative  Results 

Two  types  of  quantitative  results  were  available  from  these  tests  with 
respect  to  analyzing  unsteady  flow  behavior:  high  frequency  pressure  data 

obtained  on  the  face  of  the  flare  and  accelerometer  data  relating  the  net 
acceleration  of  the  model  and  sting  support  system  (since  the  model  is  not 
isolated) . 

Samples  of  the  four  basic  types  of  raw  pressure  and  axial  accelerometer 
data  are  shown  in  Figure  2-24.  The  meaning  of  the  exact  magnitude  of  the  axial 
accelerometer  data  is  of  little  significance  since  it  reflects  the  system 
response  of  the  entire  mode)  and  support  system;  therefore,  no  scale  is  shown. 
The  transverse  accelerations  are  difficult  to  interpret.  The  measured  approxi- 
mate maximum,  minimum  and  average  values  of  static  pressure  measured  on  the 
face  of  the  flare  are  shown  in  Table  2-4.  One  interesting  point  to  note  is 
that  the  maximum  and  minimum  values  of  pressure  for  the  metastable  type  flows 
exhibit  greater  variation  than  the  fully  pulsating  flows.  It  is  clear  from 
Figure  2-24  that  each  type  of  flow  exhibits  a unique  type  of  pressure  signal. 

The  accelerometer  signals  are  not  quite  as  unique  for  each  type  of  flow;  however, 
the  axial  accelerometer  obviously  indicates  when  strong  oscillations  are  present. 
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In  addition  to  the  raw  signal  traces  the  data  were  processed  to  obtain 
their  Power  Spectral  Densities  (PSD) . A PSD  plot  shows  the  time  averaged  ampli- 
tude squared  for  each  frequency  component  present  in  the  signal.  That  is, 
frequencies  containing  large  amounts  of  excitation  will  have  the  largest  values 
of  power  spectral  density  function.  Hence,  the  fundamental  frequency  and  its 
harmonics  axe  expected  to  exhibit  large  values  of  power  spectral  density  function. 
Figure  2-25  illustrates,  schematically,  the  nature  of  tne  pressure  and  acceler- 
ometer PSD  plots  as  related  to  the  differei t type?  of  flows  observed  visually. 

Here  again  the  pressure  PSD  plots  exhibit  characteristic  peaks  corresponding 
to  the  system  natural  frequency  shen  present,  flow  field  pulsation  frequen- 

cies . 

The  actual  pressure  a . • <,  . -'^.ometer  PSD  plots  that  were  computed  are 
shown  ir.  Figures  2-26  and  ~-2?;  not  all  runs  are  included.  All  of  these  PSD 
plots  seem  to  indicate  that  the  large  scale  flow  oscillations  occur  at  a funda- 
mental frequency  of  about  1100  cps  with  harmonics  at  about  2200  cps  and  3300  cps. 
The  fundamental  frequency  corresponds  to  a Strouhal  number,  S = vD/V,  of  about 
0.20.  There  exists  one  interesting  exception  to  the  rule,  namely,  Run  38. 

During  one  portion  of  the  run,  shown  in  Figure  2-26.38,  a strong  peak  occurred  at 


*The  pressure  and  accelerometer  oscillograms  were  generated  by  A.  Galif  of  the 
TRW  corporation. 


| 

I 2-43 


TABLE  2-4 

TABULATION  OF  HIGH  FREQUENCY  PRESSURE  DATA  (Kl) 


Run 

P 

avg 

psia 

Pmax 

psia 

Pmin 

psia 

Pto 

psia 

Tto 

•°F 

CDavg 

co 

max 

C»m)n 

1 

15 

28 

4 

260 

245 

1.69 

3.20 

.408 

2 

38 

80 

11 

575 

283 

1.94 

4.15 

.521 

4 

31 

31 

31 

607 

289 

1.49 

1.49 

1.49 

5a 

50 

117 

10 

607 

285 

2.43 

5.77 

.441 

5b 

20 

41 

15 

607 

285 

.939 

1.99 

.690 

6a 

60 

115 

15 

597 

277 

2.93 

5.67 

.703 

6b 

32 

40 

30 

597 

277 

1.56 

1.97 

1.46 

8 

39 

82 

7 

586 

265 

1.96 

4.17 

.304 

9 

41 

77 

12 

605 

284 

1.99 

3.79 

.543 

10 

44 

75 

13 

614 

282 

2.11 

3.64 

.583 

11 

40 

54 

31 

602 

275 

1.95 

2.65 

1.50 

12 

36 

40 

34 

603 

284 

1.74 

1.95 

1.65 

13 

31 

35 

26 

598 

270 

1.51 

1.71 

1.26 

14 

32 

35 

29 

599 

264 

1.56 

1.71 

1.41 

15 

28 

45 

21 

592 

275 

1.37 

2.24 

1.02 

16a 

26 

36 

21 

609 

269 

1.23 

1.73 

.985 

16b 

26 

86 

4 

609 

269 

1.23 

4.21 

.141 

17 

25 

36 

21 

607 

266 

1.19 

1.74 

.989 

18 

26 

36 

21 

600 

263 

1.25 

1.76 

1.00 

19a 

30 

48 

24 

605 

287 

1.44 

2.34 

1.14 

19b 

30 

82 

10 

605 

287 

1.44 

4.04 

.443 

20 

23 

30 

18 

600 

272 

1.10 

1.45 

.850 

21 

30 

30 

30 

606 

292 

1.44 

1.44 

1.44 

22 

32 

48 

19 

596 

277 

1.57 

2.38 

.907 

23 

35 

41 

31 

603 

266 

1.70 

2.00 

1.50 

24 

31 

45 

14 

580 

255 

1.56 

2.29 

.673 

25a(a  = 2.5°) 

22 

35 

18 

608 

288 

1.04 

1.68 

.838 

25b (a  = 5.0°) 

21 

66 

12 

608 

288 

.987 

3.22 

.540 

25c(a  = 0°) 

30 

40 

23 

600 

288 

1.43 

1.93 

1.09 

26a(a  = 0°) 

36 

41 

31 

603 

256 

1.75 

2.00 

1.50 

26b (a  = 2.5°) 

32 

45 

17 

603 

286 

1.55 

2.25 

.795 

26c(a  = 5.0°) 

28 

41 

16 

603 

286 

1.35 

2.00 

.745 

27 

24 

32 

20 

604 

302 

1.14 

1.54 

.944 

28 

18 

40 

4 

274 

288 

1.93 

4.36 

.384 

32 

34 

48 

16 

606 

276 

1.64 

2.34 

.741 

33 

17 

24 

12 

283 

268 

1.76 

2.51 

1.22 

35 

28 

39 

24 

598 

261 

1.36 

1.91 

1.16 

36 

32 

48 

19 

597 

251 

1.56 

2.37 

.905 

37 

26 

42 

20 

598 

292 

1.26 

2.07 

.954 

38 

42 

50 

34 

607 

282 

2.03 

2.43 

1.64 

Notes:  (1)  Some  records  not  available 

(2)  a and  b reflect  a change  in  character  of  the  signal 
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1)  FULLY  PUI  SATING 
EX.  RUN  2 


Figure  2-26.02.  Pressure  PSD  fnr  Run  5,  Tape  Position  01:50:30 
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Pressure  PSD  for  Run  5,  Tape  Position 
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Pressure  PSD  for  Run  8,  Tape  Position  O6:0-';30 
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Figure  2-26.09.  Pressure  PSD  for  Run  13,  Tape  Position  01:14:41 


1 ft k, a* . *?*  ' -A'i *i ill ilSVlLV^ «rt‘k5» if  a *.&Wi>.  ^ 


Frequency,  Hz 

Figure  2-26.14.  Pressure  PSD  for  Run  18,  Tape  Position  01:19:07 


PSD  PLOT 


Figure  2-26.17.  Pressure  PS?)  for^fcun^.  Tape  Position  01:17:28 


Figure  2-26.19.  Pressure  PSD  for  Run  24,  Tape  Position  01:18:04 
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Figure  2-26.21.  Pressure  PSD  for  Run  25,  Tape  Position  11:39:17 
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Figure  2-26. ?5.  Pressure  PSD  for  Run  26,  Tape  Position  11:52:02 
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Figure  2-26.30.  Pressure  PSD  for  Run  32,  Tape  Position  15:18:40.2 
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Figure  2-26.32.  Pressure  PSD  for  Run  32,  Tape  Position  15:18:44 
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Figure  2-26.40.  Pressure  PSD  for  Run  42,  Tape  Position  11:50:25 
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Figure  2-26.42.  Pressure  PSD  for  Run  42,  Tape  Position  11:50:39 
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Figure  2-26.43.  Pressure  PSD  for  Run  43,  Taps  Position  11:51:52 
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Figure  2-26.46.  Pressure  PSD  for  Run  44,  Tape  Position  12:04:51 


Pressure  PSD  for  Run  46,  Tape  Position  15:19:53 
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Figure  2-27.04.  Accelerometer  PSD  for  Run  8,  Tape  Position  06:04:30 
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Accelerometer  PSD  for  Run  9,  Tape  Position  01:12:01 
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Accelerometer  PSD  for  Run  2S,  Tape  Position 
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Figure  2-27.12.  Accelerometer  PSD  for  Run  26,  Tape  Position  11:52:12 
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Accelerometer  PSD  for  Run  42,  Tape  Position 
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Accelerometer  PSD  for  Run  43,  Tape  Position 
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Figure  2-27.22.  Accelerometer 
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about  3700  cps.  Recall  frc®  Figure  2-21,  that  this  behavior  is  apparently 
associated  with  a small  amplitude  secondary  shock  motion. 

2.6  6 Preliminary  Correlation  Results 

Two  basic  classes  of  correlations  were  attempted:  one  for  steady  sepa- 

ated  flows  and  one  for  establishing  the  limits  for  unsteady  flows.  The  correla- 
tion for  steady  separated  flows  was  intended  to  determine  the  effective  angle 
of  separation.  Figure  2-28  shows  «.  sample  of  this  type  of  calculation.  The 
observed  inclination  cf  the  tip  shock  was  measured  from  the  35  mm  photos.  From 
the  observed  shock  inclination  and  the  free  stream  Mach  number,  the  effective 
conical  angle  of  separation,  <*e££#  was  calculated  from  conical  shock  theory. 

In  all  cases  of  steady  separated  flow  the  line  determined  by  the  effective  angle 
of  separation  intercepted  the  body  at  or  near  the  shoulder  as  shown.  It  was 
also  desired  to  determine  how  the  calculated  value  of  oeff  compared  to  the  observ- 
ed edge  of  the  separated  zone  as  perceived  from  the  still  photos.  Figure  2-29 
presents  a comparison  between  the  calculated  value  of  and  the  observed  edge 

of  the  shear  layer.  In  general  the  comparison  is  close,  with  a general  tendency 
to  overestimate  the  actual  angle  of  separation  from  the  photos. 

Figure  2-30  shows  the  preliminary  attempt  to  correlate  the  data  for 
stability  limits.  The  method  of  correlation  chosen  was  similar  to  that  of 
Wood.  For  spiked  (6^  = p°)  cone  cylinders.  Wood  constructed  a correlation  for 
stability  limits  by  plotting  1/D  {or  L/D)  vs.  total  turning  angle  02.  Wood 
concluded  that  the  origin  of  the  regime  of  unsteady  flow  would  be  the  cone  angle, 

0 , that  corresponds  to  detachment  of  a conical  shock  wave.  The  data  from  the 
current  series  of  tests  is  shown  plotted  as  a Wood's  type  correlation  in  Figure 
2-30.  It  can  be  seen  that  the  correlation  for  the  limits  of  stability  seem  to 
be  fairly  distinct.  The  only  anomaly  apparent  i>  the  data  point  at  0-,  = 70°, 

1/D  = 0.5,  which  is  steady  separated  flow. 


2-121 


H 


' .w„ 


... 


TURNING  ANGLE 


u 


1 


i 

! 


SECTION  3 

ANALYTICAL  PROGRAM 

The  pulsating  flow  field  is  quite  complex.  Though  conceptually  feasible, 
an  exact  solution  of  tne  unsteady  equations  of  motion  is  not  practical  for  the 
Reynolds  numbers  of  interest.  Even  if  it  were  economically  feasible  to  solve 
the  complete  equations,  current  capability  to  model  turbulent  transport  propert- 
ies is  essentially  nonexistent  for  transient  flows  of  this  type.  Therefore, 
analytical  efforts  in  this  study  have  been  directed  toward  developing  a simpli- 
fied model  which  would  aid  in  understanding  the  cause  of  the  instability  and 
defining  the  envelope  o',  instability  onset. 

The  analysis  is  based  on  the  Korst  model  of  reattaching  separated 
turbulent  flows.  Korst* s model  has  been  used  extensively  in  analyzing  many 
separated  turbulent  boundary  layer  flows,  particularly  base  flow,  but  also 
including  the  separation  on  forward  facing  ramps  and  flares. 

Korst 's  model  is  based  upon  a reattachment  criterion,  and  it  does  not 
include  consideration  of  what  causes  the  separation.  Therefore,  it  must  be 
combined  with  a model  of  the  separation  region  before  the  separated/reattaching 
flow  field  problem  can  be  solved  in  toto.  Compare,  for  a moment,  separation 
over  a rearward  facing  step  with  that  over  a forward  facing  ramp  (Figure  3-1) . 

As  noted  above,  Korst's  model  relates  the  flow  in  rer  >ns  2 and  3,  using  simpli- 
fying assumptions  about  the  details  of  the  reattachment  process  (R  in  Figure  3-1) . 
For  the  rearward  facing  step,  the  flow  external  to  the  shear  layer  in  region 
2 and  the  reattachment  point  R are  unique  functions  of  the  approach  Mach  number 
and  turning  angle  at  the  corner.  The  separation  length  is  uniquely  determined 
by  the  angular  geometry  set  by  the  separation/reattachment  process  and  the 
step  height,  H.  For  the  forward  facing  ramp,  if  there  is  sufficient  running 
length  ahead  of  the  corner,  the  edge  conditions  in  region  2 are  again  uniquely 
determined  by  the  flow  in  region  1 and  the  shock  induced  turning  angle  between 
regions  1 and  2.  In  such  a case  the  turning  from  region  1 to  2 is  governed  by 
a free  interaction  process.*  There  is  no  corner  imposing  the  separation  point, 
and  there  is  no  length  scale  in  Korst's  analysis  with  which  to  select  the 


See  below  for  a discussion  of  free  vs.  forced  interaction. 
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a)  Rearward  Facing  Step 


b)  Forward  Facing  Ranp 


Figure  3-1.  Schematic  of  Separation  over  a Rearward 
Facing  Step  and  a Forward  Facing  Ramp 

separation  length.  I - . •<  '"v,ns.  it  is  necessary  to  include 

- Co*icidtiou  ror  separation  length  to  solve  fcne  av.’  faction  problem  of 

the  forward  facing  ramp. 

There  is  another  aspect  in  which  the  forward  facing  ramp  differs  from 
the  rearward  facing  step.  In  the  latter  case  the  flow  definitely  separates  at 
the  corner  and  reattaches  on  the  recompression  surface.  But  for  the  forward 
facing  ramp,  the  flow  may  not  separate  at  all,  or  it  may  not  be  fully  separated. 
Thus,  the  Korst  type  analysis  is  applicable  when  the  flow  is  fully  separated 
in  the  sense  that  the  velocity  profile  in  the  shear  layer  has  reached  a develop- 
ed shape  and  remains  unchanged  in  shape,  that  is,  the  velocity  profiles  are 
similar. 

As  with  many  flow  models  which  are  relatively  simple  physically,  details 
of  Korst  type  analyses  of  reattaching  flow  are  somewhat  complex.  Details  of 
the  present  analyses  are  included  in  the  Appendix. 

To  this  point  the  discussion  has  been  about  steady  flow,  and  a Korst 
type  analysis  in  its  basic  form  provides  no  information  about  the  transient 
in  an  unsteady  flo*-  process.  On  the  other  hand,  one  may  anticipate  that  such 
an  analysis  may  approximately  model  physical  interactions  which  define  boundaries 
between  different  flow  regimes.  That  is,  the  Korst  type  analyses  are,  in  a 
sense,  analyses  of  limiting  cases  which  may  provide  insight  into  the  nature  of 
the  boundaries  between  different  flow  regimes. 


We  now  set  aside  the  rearward  facing  step  problem  and  concentrate  on 
the  question  of  the  forward  facing  ramp.  With  regard  to  separation,  the 
question  arises  as  to  what  occurs  when  the  body  surface  extends  upstream  a 
distance  less  than  the  length,  1 , defined  by  the  free  interaction  geometry. 
This  corresponds  to  the  condition  where  the  separation  occurs  at  or  near  the 
leading  edge  of  the  body.  In  such  a case,  common  plateau  pressure  criteria 
for  establishing  the  angle  of  separation,  <J>,  are  inapplicable  since  there  is 
not  an  established  boundary  layer  approaching  the  point  of  separation.  This, 
therefore,  indicates  that  no  longer  can  the  angular  geometry  at  the  separation 
point  be  set  by  the  local  conditions  at  the  point  of  separation,  this  being 
the  definition  of  a free  interaction,  but  instead  is  set  by  a downstream  re- 
attachment criteria.  The  latter  type  of  separation  is  termed  a forced  inter- 
action as  opposed  to  a free  interaction. 

Recall  that  the  Korst  reattachmer.t  model  provides  no  information  about 
the  separation  process.  For  free  interaction  at  separation,  the  plateau  pres- 
sure coefficient  for  region  2 can  be  correlated  approximately  with  the  relation 
ship 


8 

1 + 10  Mx 


(3-1) 


This  correlation  reflects  the  fact  that  for  fully  turbulent  separations,  the 
plateau  pressure  rise  is  largely  independent  of  the  Reynolds  number.  Thus, 
for  given  M^,  the  pressure  in  region  2 can  be  obtained  from  the  Equation  (3-1) . 
The  separation  angle  0 and  all  edge  conditions  in  region  2 can  then  be  calculat 
ed  from  oblique  shock  relatione.  Then  the  Korst  analysis  can  be  applied  to 
determine  the  unique  ramp  angle  which  satisfies  that  reattachment  criterion. 

To  illustrate  some  of  these  points  further,  the  sequence  of  events 
defining  the  transition  from  a free  interaction  to  a forced  interaction  is 
schematized  in  Figure  3-2.  Figure  3-2a  illustrates  a Type  I free  interaction 
in  which  the  separation  has  not  reached  its  full  extent.  In  this  case  the 
ramp  angle,  8,  is  less  than  the  critical  angle  for  reattachment  computed  from 
Korst 's  theory.  That  is,  8 represents  the  ramp  angle  associated  with  a 
large  separated  region  rn  which  the  shear  layer  has  become  fulxy  developed. 

(If  the  separated  region  is  very  small,  then  free  interaction  may  not  exist 
because  the  separation  and  reattachment  processes  may  directly  affect  each 
other.)  For  the  Type  I free  interaction,  there  is  sufficient  body  length  so 
that  a fully  developed  boundary  layer  is  established  ahead  of  the  point  of 
separation,  and  the  angular  geometry  at  the  separation  point  is  determined 
independent  of  the  reattachment  process,  and,  for  example,  a plateau  pressure 
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Figure  3-2. 
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correlation  such  as  that  represented  by  Equation  3-1  is  appropriate.  The  length 
of  the  separated  region  is  governed  by  the  height  of  the  approaching  boundary 
layer  and  the  ramp  angle.  Because  the  shear  layer  has  not  become  fully  develop- 
ed, the  Korst  theory  of  reattachment  is  not  applicable. 

Figure  3-2b  illustrates  the  second  type  of  free  interaction  separation. 

In  this  instance  the  ranp  angle  is  increased  to  a value  equal  to  that  which 
yields  a fully  separated  region  having  a fully  developed  shear  layer.  The  Korst 
reattachment  model  is  applicable;  hence,  the  critical  ramp  angle,  8cri^,  corres- 
ponds tc  th“  5nimum  angle  for  which  the  shear  layer  is  fully  developed.  Again 
since  an  established  boundary  layer  preceeds  the  point  of  separation,  the  geo- 
metry at  the  point  of  separation  is  governed  primarily  by  the  approach  Mach 
number  and  the  height  of  the  approaching  boundary  layer.  The  turning  at  re- 
attachment  is  governed  by  the  Korst  reattachment  model  and  oblique  shock  relations. 
This  type  oC  flat-;  separation/reattachment  can  be  fully  analyzed  within  the  con- 
text of  the  theory  developed  herein. 

Now  suppose  that  the  ramp  angle  is  increased  further  to  some  angle,  B, 
greater  than  °-cr-i ?s  shown  in  Figure  3-2c.  This  represents  the  third  type 
of  free  interaction.  In  f-his  case,  since  the  shear  layer  has  already  reached  its 
fully  develop-^  shape,  and  since  the  ramp  angle  is  greater  than  that  governed  by 
the  Kor3t  re attachment  model,  reattachment  on  the  ramp  face  is  not  possible. 

Hence,  the  separation  point  is  forced  to  move  upstream  so  that  reattachment 
occurs  at  the  shoulder  of  the  ramp.  Since  the  approach  Mach  number  dominates 
the  angular  geometry  at  separation  (e.f.  the  plateau  pressure  coefficient 
correlation,  Equatior  _ 1) , the  length  of  the  separation  zone  is  no  longer 
dependent  upon  the  height  of  the  approaching  boundary  layer  but,  instead,  upon 
the  height  of  the  ramp.  Zukoski7  shows  this  to  be  the  case  for  various  forward 
facing  steps;  that  is,  the  flow  geometry  now  behaves  as  though  the  separation 
occurred  ov'-r  a forward  facing  step. 


f The  three  situations  depicted  in  Figures  3-2a-c  have  one  thing  in  common, 

r the  incident  boundary  layer  develops  over  a long  length  and  there  are  no  geo- 

5 metrical  considerations  restricting  the  separation  point.  That  is,  the  separ- 

r 

5 ation  location  and  geometry  are  free  to  adjust  themselves  according  to  local 

i flow  field  interactions.  Obviously,  separation  could  be  forced  to  occur  at 

J some  point  s downstream  of  a free  interaction  separation  point  simply  by  removing 

all  of  the  wetted  surface  upstream  of  s.  Two  such  forced  separation  situations 
? are  of  interest  in  this  study,  differing  in  the  nature  of  the  flow  at  reattach- 

> ment. 
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In  the  first  such  forced  separation  case,  the  separated  shear  layer 
reattaches  on  the  ramp  surface  upstream  of  the  shoulder.  Figure  3-2d.  There- 
fore, Korst’s  model  is  applicable  to  the  reattachment  process,  and  the  theory 
for  this  type  of  forced  interaction  predicts  the  angle  of  separation  in  addition 
to  the  total  amount  of  turning  possible.  Hence,  there  is  a continuous  variation 
of  Type  I forced  interactions  possible  up  to  a maximum  value  of  ramp  angle, 

3^  ^ . Also  as  the  flare  angle,  3,  decreases,  there  is  an  obvious  lower 

limitm!^r  the  Type  I forced  interaction  result,  that  being  when  3 * 3 j_t, 
the  free  interaction  value  corresponding  to  the  Korst  limiting  case.  When 
3 < 3cr£t  the  flow  no  longer  separates  from  the  leading  edge  but  returns  to 
one  or  the  free  interaction  models  shown  in  Figures  3-2a,  b. 


Jf 

I 

$ 


i 


I ! 
E ! 

I • 


L- 


In  the  second  type  of  forced  interaction,  the  reattachment  point  is 
also  determined  by  geometrical  conditions  rather  than  local  flow  field  inter- 
actions. We  designate  S^rit  as  t*ie  va3-ue  of  3 for  which  the  separated 

shear  layer  reattaches  at  th^Shoulder , Figure  3-2e.  Certainly  an  upper  bound 
on  f^rit  the  ramp  angle  corresponding  to  sonic  flow  behind  the  reattach- 

ment sho8*?  The  separated  shear  layer  geometry  for  the  Type  II  forced  inter- 
action is  set  by  an  equivalent  ramp  drawn  from  the  leading  edge  to  the  aft 
shoulder. 

In  order  to  validate  the  hypothesis  that  there  exists  a limiting  ramp 
angle  for  a free  interaction  reattachment,  a study  was  initiated  utilizing 
existing  separation  data.  A detailed  investigation  of  this  study  is  presented 
in  the  Appendix.  In  summary,  the  ratio  of  separation  length  to  approach  bound- 
ary layer  height  was  computed  as  a function  of  ramp  angle.  The  experimental 
data  is  generally  presented  in  the  form  of  static  pressure  ratio  distributions. 

By  taking  the  plateau  pressure  ratio  and  assuming  that  the  separated  region 
acts  essentially  as  an  equivalent  ramp  surface,  a separation  angle  was  computed 
using  oblique  shock  relations.  Then  by  measuring  the  distance,  13,  from  the 
ramp  corner  to  the  point  of  separation  and  applying  appropriate  trigonometric 
relationships,  the  length  of  the  separated  shear  layer,  lgep,  Figure  3-1,  was 
calculated.  According  to  the  preceeding  discussion  for  a free  interaction  re- 
attachment on  the  ramp  face.  Type  II,  there  should  exist  a limiting  ramp  angle, 
Bcrifc,  such  that  the  ratio  igep/^0  approaches  this  limit.  Figure  3-3  illustrates 
the  expected  behavior  of  the  data.  Actual  data  are  not  as  well  behaved  as 
illustrated  in  Figure  3-3.  However,  examination  of  the  results  shown  in  the 
Appendix  illustrates  that  the  general  concept  shown  in  Figure  3-3  is  likely  to 
be  valid.  Naturally,  owing  to  the  indirect  means  in  which  the  separation  data 
were  obtained  from  the  pressure  distributions,  considering  that  the  approaching 
boundary  layer  is  affected  by  the  presence  of  boundary  layer  trips,  and,  in 


3-6 


-j  • ‘g».-wa»  <rcr,-  ;-.<•— r:  ~r-~—- 


'"incip  'crit 


Figure  3-3.  Schematic  of  Hypothesized  limit  of 
Separation  Length 

addition,  allowing  for  a certain  amount  o*  inexactness  in  the  data,  the  scatter 
in  the  results  obtained  is  not  unexpected. 

The  preceading  discussion  pertaining  to  the  character  of  separated  flows 
applies  to  the  hypothesized  theoretical  model  proposed  herein.  Portions  of  this 
theory  are  taken  from  observed  experimental  evidence.  This  is  not  to  say  that 
the  theory  is  all  encompassing.  For  instance  the  theory  cannot  predict  when 
separation  occurs;  on  the  other  hand,  the  question  may  be  asked  as  to  whether 
the  theory  can  be  used  to  predict  the  onset  of  flow  instabilities  for  separated 
flows.  For  the  case  of  free  interaction  type  separations,  experimental  evidence 
indicates  that,  if  6cr^t  is  exceeded  at  a given  free  stream  condition,  the  flow 
may  become  unstable. 

Figure  3-4  illustrates  this  point  more  fully  by  comparing  the  present 

theory  for  predicting  6 . for  a two-dimensional  free  interaction  with  experi- 

(5)  (6) 

mental  data  gathered  by  Kuehr.'  'et  al  . This  data  indicate  that.,  for  a free 
interaction  type  separation,  some  type  of  flow  anomaly  is  likely  to  result  if 
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Figure  3-4.  Theory  Compared  with  Observed  Flow  Behavior  on 
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the  ramp  angle  exceeds  $cr^t  as  predicted  by  Korst's  theory.  These  anomalies, 
as  observed  by  Kuehn,  et.  al.,  may  take  one  of  three  forms: 

i)  The  flow  may  be  grossly  unsteady  with  large  scale  oscillations 

ii)  If  the  flow  is  not  separated,  a hysterisis  in  the  pressure 
distribution  may  be  exhibited  with  increasing  and  decreasing 
Mach  number  variation,  and 

iii)  If  the  ramp  angle  exceeds  the  critical  value,  &cr^t*  the  flow 
may  jump  from  an  attached  flow  to  a fully  separated  flow. 

These  behavior  patterns  indicate  that  the  boundary  defined  by  is  signi- 

ficant in  predicting  the  nature  of  the  stability  of  the  flow  separation. 

The  question  now  arises  as  to  how  the  Korst  theory  can  be  used  to  gain 
insight  into  the  nature  of  flow  separations  occurring  on  the  nosetip  shapes  of 
interest  and  how  this  information  can  be  related  to  the  prediction  of  the  limits 
of  flow  instability  on  these  shapes.  For  shapes  of  interest,  if  02  falls  within 
the  limits  of  allowable  reattached  solutions , then  the  method  can  be  used  to 
predict  the  free  or  forced  interactions  which  exist,  and  the  flow  should  be 
steady  (Figures  3-2b,  c,  d) . The  results  of  the  analysis  for  forced  interaction 
(Figure  A6-4  of  the  Appendix)  indicate  that  the  values  of  02  tested  in  this 
program  exceed  the  value  of  so  that  the  reattachment  should  be  forced 

to  the  shoulder.  Also,  the  tip  lengths  tested  art  considerably  less  than  lengths 
consistent  with  free  interaction  at  separation  for  such  large  values  of  02; 
therefore,  it  is  concluded  than  in  these  tests  the  separated  flows  should  be 
Type  II  forced  interaction  flows.  Figure  3-2e.  In  fact  this  was  observed  to 
be  the  case  for  all  steady  separated  flows  in  the  test  program. 

The  next  question  that  is  posed  is  how  the  Korst  model  can  be  used  to 
predict  the  limits  of  flow  stability,  if  at  all.  Consider  the  shapes  of  interest 
shown  in  Figure  3-5.  Figure  3-5a  depicts  the  case  discussed  in  the  previous 
paragraph  relating  to  the  steady  forced  interaction  which  results  in  the  flow 
reattaching  at  or  near  the  shoulder.  However,  if  the  geometry  is  such  that 
the  tip  protrudes  through  the  bow  shock  that  would  exist  if  the  tip  were  not 
present,  and  if,  in  addition,  the  ratio  of  1/D  as,  shown  in  the  sketch  is  such 
that  an  attached  conical  oblique  shock  could  not  be  r<  alized,  then  it  is  ex- 
pected that  the  flow  will  be  unstable,  Figure  3 -5b.  That  is  to  say,  if  the 
effective  conical  body  angle,  created  by  the  separation  exceeds  the 

maximum  value  for  the  occurence  of  an  attached  oblique  shock  then  flow  instabil- 
ity will  occur.  The  cause  is  illustrated  in  Figure  3-5c  and  3-5d.  If  the 
separation  were  to  occur  from  the  tip  to  the  shoulder,  a strong  shock  would 
have  to  appear  as  shown  in  Figure  3-5c  (assuming  the  flow  is  steady) . However, 
this  steady  state  can  not  persist  because  the  large  separation  bubble  required 
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cannot  be  sustained  by  the  subsonic  flow  behind  the  strong  shock.  The  bubble 
must  collapse,  as  shown  in  Figure  3-5d.  However  the  flow/geometry  conditions 
are  incompatible  with  the  limited  separation  region  (Figure  3-5d) , and  the 
separation  point  moves  back  toward  the  tip,  leading  to  the  observed  oscillatory 
flow.  This  does  not  preclude,  however,  the  existence  of  steady  flow  with  a 
completely  attached  boundary  layer  and  a shock  structure  similar  to  that  in 
Figure  3-5d. 
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Figure  3-5.  Limits  of  Instability  as  Predicted  from  Theory 

In  the  current  experiments  it  was  shown  that  when  the  flow  is  separated 
and  oef£  exceeds  the  maximum  value  for  an  attached  shock,  flow  instabilities  did 
occur.  This  seems  to  be  the  extent  to  which  the  current  theory  can  be  used 
to  predict  the  occurence  of  flow  instability.  That  is,  it  was  observed  in  the 
current  test  series  that  a range  of  strong  instabilities  exists  for  instances 
where  the  1/D  did  not  indicate  that  c_m  „ would  be  exceeded.  The  cause  of  tnis 
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SECTION  4 

STABILITY  CORRELATION  CONSIDERATIONS 
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A comparison  of  the  current  correlation  technique  was  made  with  the  data 
taken  by  Wood  at  Mach  10  as  shown  in  Figure  4-1.  It  is  expected  that  above 
M =*  5 there  will  not  exist  a strong  Mach  number  dependency  upon  stability  limits. 
This  is  based  upon  the  fact  that  the  conical  shock  detachment  angle  does  not 

g 

vary  greatly  with  Mach  number  above  M = 5.  Secondly,  as  previously  mentioned, 
the  upper  limit  of  the  stability  regime  is  associated  with  separated  flow  attach- 
ing at  the  shoulder,  and  for  M ^ 5 the  shock  structure  is  relatively  insensitive 
to  changes  in  Mach  number.  The  comparison  of  the  regimes  of  unsteady  flow  for 
the  current,  tests  and  Wood's  data  illustrates  a drastic  difference  in  the  extent 
of  the  limits.  Figure  4-1.  In  addition  the  data  point  shown  for  an  AVCO/VKI  test 
indicates  that  according  to  Wood's  data  the  flow  should  have  been  unsteady,  but 
in  fact  the  flow  was  steady.*  The  Series  E data  ooints,  which  were  both  unsteady 
are  also  shown  in  Figure  4-1. 

The  one  obvious  difference  between  Wood's  models  and  the  present  series 
qc  models  is  the  difference  in  the  tip  angle;  that  is.  Wood  only  tested  models 
with  spike-like  protrusions,  0^  = 0°,  and  in  the  present  series  of  tests  the 
tip  angle  varied  from  20°  to  40°.  It  seems  clear  that  the  tip  angle  probably 
has  a strong  influence  upon  the  stability  limits.  One  reason  for  this  postula- 
tion is  that,  as  mentioned  previously,  as  the  tip  angle  increases  and  approaches 
the  line  tangent  from  the  nose  to  the  shoulder  then  the  body  becomes  a plain 
cone  cylinder  and  is  obviously  steady.  At  the  other  extreme  when  0^  is  small, 
that  is,  when  the  tip  approaches  a spike  configuration,  then  for  a forced  separ- 
ation the  separation  bubble  volume  is  large.  Hence,  a stability  condition  can 
be  hypothesized  as  follows.  If  a forced  separation  exists  and  the  bubble  is 
large,  assume  that  the  dividing  streamline  is  perturbed  from  being  straight  and 
instead  assumes  a disturbed  sinusoidal  condition.  The  question  arises  as  to 
whether  or  not  this  disturbance  will  amplify  and  cause  a flow  instability  to 
exist.  In  a very  general  sense  the  forces  which  will  tend  to  nullify  the  dis- 
turbance are  the  momentum  in  the  shear  layer  and  the  restoring  force  of  the 
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bubble.  Viscous  effects  in  the  shear  layer  will  tend  to  retard  shear  layer 
momentum.  The  restoring  force  of  the  bubble  is  inversly  proportional  to  the 
size  of  the  bubble;  hence,  one  would  expect  that  for  small  9^  the  flow  is  more 
likely  to  become  unsteady.  This  is  a broad  brush  treatment  of  a very  difficult 
problem;  however,  the  basic  precepts  seem  very  reasonable. 

It  was,  therefore,  concluded  that  6,  must  play  a significant  role  in 

creating  the  potential  for  unsteady  flow.  The  body  of  data  from  the  current 

test  series  was  replotted  along  with  Wood's  data  in  a new  correlation  format 

in  Figure  4-2.  Here  another  dimension  has  been  added  to  the  correlation,  namely, 

planes  of  constant  6^.  In  place  of  02  used  by  Wood,  it  seems  appropriate  to  use 

A 6 = e2  - 61  which  leaves  Wood's  data  unchanged  since  0^  =0°.  In  addition 

the  left  hand  boundary  which  was  previously  based  upon  the  critical  angle  for 

shock  detachment,  6_  , for  Wood's  data  was  extended  to  correspond  to  the  max- 

^max 

imum  combination  of  angles  0 =1(9.,  + & 0)]  to  yield  shock  detachment  via 

zndx  X max 

a conical  shock  plus  a two  dimensional  shock.  This  left  hand  limit  was  establish- 
ed for  M = 5 but  should  not  change  significantly  for  higher  Mach  numbers. 

The  result  of  this  correlation  is  a volume  of  instability;  that  is,  any 
points  falling  inside  the  volume  should  be  unstable.  One  of  the  particularly 
satisfying  aspects  of  this  correlation  is  that  the  previously  mentioned  AVCO/VKI 
data  point  falls  outside  of  the  volume  of  instability  and  explains  why  the  flow 
was  steady  in  that  test.  The  volume  shown  in  Figure  4-2  also  makes  it  apparent 
that  there  is  a significant  increase  in  the  limits  of  the  unstable  flow  regime 
as  0 ^ decreases. 


4-3 


MODIFIED  WOOD(m.»io) 


SECTION  5 


CONCLUDING  REMARKS 


The  fluid  mechanical  interactions  associated  with  the  flow  fields  of 
interest  in  this  study  are  quite  complex.  It  is  certainly  within  existing 
technical  capability  to  develop  computer  codes  to  obtain  exact  solutions  of 
the  full  Navier-Stokes  equations  of  motion  (though  the  problems  associated 
with  accurately  representing  the  transient  shock  structure  are  formidable) . 

Such  an  approach  would  enable  us  to  study  the  detailed  interactions  which 
lead  to  instability  onset  and  to  parametrically  vary  key  shape  and  flight 
configurations  to  determine  the  dependence  of  instability  onset.  Each  such 
calculation  would,  in  itself,  require  an  exhorbitant  amount  of  computer  time, 
making  such  an  approach  economically  unacceptable . In  addition,  we  have  very 
little  understanding  of  how  to  model  turbulent  transport  properties  for  tran- 
sient flow,  and  the  regimes  of  interest  in  flight  are  for  transitional  and 
fully  turbulent  boundary  layers.  Therefore,  this  program  is  oriented  toward 
a semi-empirical  approach  to  solve  the  problem  of  pulsating  flow  on  ablated 
shapes. 

The  result  is  a collection  of  experimental  data  and  empircial,  semi- 
empirical  and  approximate  analytical  models  which  are  patched  together  to 
form  a "solution"  to  the  problem.  As  is  often  the  case,  the  simpler,  engi- 
neering approach  is  more  piecemeal  them  an  exact  numerical  solution  would  be, 
because  the  engineering  solution  is  a synthesis  of  many  fragments.  In  such  a 
situation  it  is  easy  to  focus  on  the  parts  and  to  lose  sight  of  the  whole  of 
the  model.  Therefore,  we  conclude  with  a synopsis  of  the  major  conclusions 
of  this  study. 

The  overall  objective  is  to  define  the  envelope  of  instability  onset 
for  ablated  shapes  of  interest  in  flight  environments.  There  is  no  data  to 
define  shapes  associated  with  flight  !to  the  required  precision  for  our  needs. 
The  flow  field/model  shape  interaction  is  qualitatively  the  same  fcr  wind 
tunnel  conditions  as  it  is  for  flight,  and  Mach  number  effects  are  relatively 
weak  and  can  be  adequately  accounted  tor  with  approximate  analyses.  Therefore, 
the  solution  is  based  cn  ground  test  data  of  low  temperature  ablator  shape 
change  experiments. 


The  highlights  of  the  study  are: 

• Realistic  ablated  shapes  are  represented  by  a simplified  but 
adequate  model  which  has  been  described  by  two  angles,  three 
radii,  and  a length. 

• The  three  parameters  most  important  to  instability  onset  are 

the  two  angles,  0^  and  02»  and  the  length  to  diameter  ratio,  i/D. 

• The  flare  and  shoulder  radii,  and  Rg,  are  of  secondary  impor- 
tance for  the  range  of  shapes  likely  to  be  encountered  in  realis- 
tic trajectories. 

• The  nose  radius  appears  to  be  of  intermediate,  but  closer  to 
secondary  than  primary,  importance.  The  phenomenology  of  the 
nose  radius  effect  is  unclear. 

• As  long  as  the  flow  is  fully  turbulent,  Reynolds  number  varia- 
tions have  a weak  effect. 

• For  H i 5,  Mach  number  effects  should  be  relatively  small  and 
predictable. 

• For  shapes  and  conditions  of  interest,  separation  of  the  boundary 
layer  (when  it  separates)  is  of  the  forced  type  and  the  shear 
layer  reattaches  at  or  very  near  the  shoulder.  Therefore,  free 
interaction  theory  results  cannot  be  used  directly  to  predict 
the  separated  flow. 

• As  noted  by  Maull  and  Wood,  the  inviscid,  conical  detachment 
angle  for  the  recompression  cone  (02)  is  important  in  determining 
the  envelope  of  instability  onset. 

• Wood's  envelope  of  instability  onset  was  modified  and  extended  to 
realistic  ablated  nosetip  shapes.  The  forecone  angle,  6^,  is  a 
strong  parameter  in  this  correlation,  and  as  6^  increases,  the 
range  of  ©2  and  1/D  for  which  the  flow  field  pulsates  decreases 
rapidly. 
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APPENDIX  A 

ANALYSIS  OF  STABILITY  OF  FLOW  SEPARATIONS  ON 
EICONIC  SURFACES 


Larry  Cooper 
Thomas  J.  Dahm 


LIST  OF  SYMBOLS  FOR  APPENDIX  A 


Crocco  number  at  edge  of  shear  layer 


pressure  rise  coefficient 


error  function 


dimensionless  bleed  rate 


I,*1? 


integral  properties  of  the  shear  layer 


l ,i 
sep 


length  of  shear  layer 


Mach  number 


bleed  flow  per  unit  width 


static  pressure 


total  pressure 


static  temperature 


total  temperature 


velocity  in  direction  of  reference  coordinates 


velocity  transverse  to  reference  coordinates 


intrinsic  coordinates 


reference  coordinates 


angle  between  reference  axis  and  reattac'nment  streamline 


LIST  OF  SYMBOLS  FOR  APPENDIX  A (Continued) 


angle  separation  makea  with  flow  in  Region  1 


ramp  or  total  turning  angle 


ramp  angle  for  fully  developed  shear  layer 


INCIP 


ramp  angle  for  which  separation  just  begins 


maximum  ramp  angle  for  forced  separation 


height  of  approaching  boundary  layer 


inviscid  flow  turn  from  1 to  2 


inviscid  flow  turn  from  2 to  3 


distance  from  corner  to  point  of  separation 


equal  to  ym  in  magnitude 


dimensionless  coordinate 


turbulent  mixing  parameter 


density 


eddy  viscosity 


angle  of  tip  surface 


shock  angle  off  tip 


shock  angle  off  flare  or  ramp 


ratio  of  specific  heats 


static  temperature  ratio,  T/T2 
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SECTION  A1 


INTRODUCTION 

The  possibility  of  encountering  unsteady  flow  oscillations  about  biconic 
and  spiked  bodies  has  been  noted  in  the  literature.  Several  experimenters  such 
as  Wood^A3^,  Maull^A^,  and  Mair^A3^  observed  these  oscillations  on  wind  tunnel 
models  and,  to  a limited  extent,  correlated  the  bounds  of  steady  and  oscillatory 
flow  based  upon  the  geometry  of  the  body.  Few  attempts  were  made  to  model  the 
mechanism  which  triggers  the  unsteady  condition.  Maull  suggested  that  the 
instabilities  are  triggered  by  separations  which  lead  to  strong  family  reattach- 
ment shocks.  He  suggested  that  this  would  lead  to  a growth  and  collapse  of  the 
separated  zone.  The  present  paper  offers  a proposed  model  which  is  intended 
to  predict  the  separation  characteristics  about  a biconic  shape  and  indicate 
when  oscillatory  flow  is  likely  to  occur. 
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SECTION  A2 
APPROACH 


The  approach  to  the  problem  of  studying  the  onset  of  flow  ins tabili ties 
was  to  first  look  at  the  possible  causes  for  unsteadiness  rather  than  to  attempt 
a full  model  of  the  oscillatory  process  itself.  Frort  a design  standpoint  the 
critical  question  is  when  the  oscillations  occur  and  how  they  can  be  avoided. 


At  present  the  general  consensus  is  that  the  flow  oscillations  are  re- 


lated to  changes  in  the  effective  geometry  of  the  body  due  to  growth  of  the 
(an  (a2)  (A3) 

separation  bubble.  The  bubble  is  thought  to  grow  to  some  critical 


state  at  which  time  further  growth  cannot  be  sustained  and  the  bubble  suffers 
a collapse  due  to  outflow  near  the  shoulder.  The  exact  mechanism  of  bubble 
growth  and  collapse  is  a matter  of  speculation  at  tnis  point. 


The  approach  taken  in  the  present  study  was  to  examine  the  separation 
and  reattachment  mechanisms  for  the  steady  flow  case  and  to  observe  when  these 
steady  flow  mechanisms  become  violated.  It  is  hypothesized  that  when  these 
critical  steady  separated  conditions  sure  violated  that  a readjustment  in  the 
flow  relating  to  either  steady  or  quasi-pulsating  reattachment  at  the  shoulder 
will  occur  or  the  flow  wiil  become  completely  oscillatory. 


The  present  flow  model  consists  of  a mechanisms  for  flow  separation, 
either  free- interaction  or  forced  sepsuration,  a free  shear  flow  analysis  of  the 
sepsurated  zone  based  upon  the  theory  of  Ko:  st,  and  reattachment  turning  condi- 
tions relating  to  oblique  shock  theory. 
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SECTION  A3 


ii 


A3. 1 BOUNDARY  LAYER  SEPARATION 

There  are  several  references  in  the  literature  which  discuss,  in  detail, 
what  is  known  about  the  mechanism  of  turbulent  flow  separation  induced  by  steps 
and  rai.p3  in  supersonic  flow.  Suffice  it  to  say  that  the  separation  is  caused 
by  the  fact  that  the  momentum  of  the  boundary  layer  flow  cannot  negotiate  a 
rapid  rise  in  pressure  associated  with  the  compression  at  a corner. 

figure  A3-1  illustrates  the  salient  features  of  the  type  of  flow 

separation  considered.  At  a given  value  of  approach  Mach  number,  K^,  there  will 

exist  some  critical  ramp  angle,  81NCIP»  which  corresponds  to  the  appearance  of 

an  inflection  point  in  the  static  pressure  rise  as  measured  at  the  wall  surface 

(figure  A3- la) . As  the  angle  8 is  increased  to  values  greater  than  3INCIP» 

holding  fixed,  the  extent  of  the  separated  zone  continues  to  grow.  Figure 

A?-lb  illustrates  a typical  pressure  distribution  for  a fully  separated  region. 

The  pressure  rises  from  its  initial  value  up  to  a so-called  plateau  level. 

Experiments  have  shown  that  the  rise  to  the  plateau  level  is  mainly  Mach  number 

dependent  and  that  the  plateau  level  agree.,  well  with  the  inviscid  pressure 

in  Region  2.  During  this  phase  of  separation  the  length  of  the  separated  region, 

*-s„p/  correlates  well  with  the  height  of  the  approaching  boundary  layer,  6q, 

and  the  plateau  pressure  coefficient,  Cp^  ^ • However,  it  is  well  known  that, 

foi  a forward  facing  step  with  sufficient  running  length  ahead  of  the  step, 

(A4) 

the  separated  flow  will  always  attach  at  the  shoulder  . In  this  case  the 
length  of  tne  separated  reqion  no  longer  is  dependent  upen  the  height  of  the 
approaching  boundary  layer  but  instead  upon  the  height  or  the  and  the 

plateau  pressure  coefficient.  The  question  then  arises  whether  there  exists 
sc«oe  critical  ramp  angle,  6 ^t,  Figure  A3  2,  that  corresponds  to  an  asymptotic 
limit  for  . It  is  supposed  that  at  this  value  of  8 the  separated  shear 

layer  reaches  a similar  profile  and  achieves  maximum  allowable  reattachment 
pressure  on  the  face  of  the  r_mp.  A further  increase  in  8 will  then  drive  the 
reattachment  toward  the  shoulder  while  the  shear  layer  separation  ancle  and  the 
plateau  pressure  coefficient  are  unchanged.  Existing  data  have  been  analyzed 
to  verify  this  hypothesis.  The  results  are  summarized  in  Section  A3, 2- 
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Figure  A3-2.  Schematic  of  Hypothesized  Limit 
of  Separation  Length 
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EXTENT  OF  SEPARATION  LEGNTHS  FOR  COMPRESSION  CORNERS 


The  literature  contains  a limited  number  of  attempts  to  correlate 
separation  lengths.  None  of  these  attempts  to  correlate  ^SEP  is  wholly  extensive 
or  successful.  The  available  data  consists  mainly  of  two-dimensional  and  nearly 
two-dimensional  (flared  surfaces  far  frcm  the  axis)  pressure  distributions. 
Therefore,  little  or  no  data  is  available  which  considers  a strong  radial  effect. 
In  addition  the  literature  is  nearly  vcid  of  consideration  of  flow  separations 
such  tiiat  the  requi red  turning  at  reattachment  exceeds  the  critical  turning 
angle  for  an  attached  oblique  shock.  The  small  amount  of  data  that  is  available 
indicates  that  the  flow  may  be  unstable  if  this  condition  is  exceeded.  A third 
type  of  data  that  is  lacking  is  that  which  considers  separations  at  the  leading 
edge  in  which  the  free  interaction  would  no  longei  allow  reattachment  at  the 
ramp  shoulder;  this  would  be  termed  a forced  separation. 

To  correlate  turbulent  separation  lengths  at  supersonic  Mach  numbers 
for  ramp-induced  separation,  we  define  the  point  of  separation  as  shown  in 
Figure  A3- lb.  If  the  plateau  pressure  level  is  not  clearly  definable,  an 
empirical  plateau  coefficient  such  as 


3.  2 

Cp  = g-—  — r”ip  (Reference  A5) 

is  used  to  define  the  piateau  level.  Given  the  plateau  level  either  by  inspec- 
tion of  the  experimental  pressure  distribution  or  from  the  above  expression , 
the  angle  of  separatic  . car.  be  determined  from  the  two-dimensional  oblique  shock 
relations . 

6 12  = 5 12  jMl' 


With  the  angle  of  separation,  and  the  measured  value  of  AX,  the  separation 

length,  J-SEP*  can  be  determined  from  geometric  considerations. 


‘'SEP  £SEP  jAX'  aSEP'  6 1 


This  procedure  defines  a uniform  method  of  reducing  the  pressure  distribution, 
data  into  the  desired  form  The  data  sources.  References  A6-A14,  are  shown  in 
Figure  A3. 3.  The  data  is  shown  plotted  in  Figures  A3-4  to  A3-8.  Also  shown 
is  the  two-dimensional  free-interaction  critical  angle  for  reattachment.  That 
is,  the  ramp  angle  that  produces  a separation  which  extends  sufficiently  so 
that  the  velocity  profile  in  the  shear  layer  has  become  similar  (see  Section  3.0 
in  body  of  report).  . . 
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SECTION  A4 


TURBULENT  MIXING  ANALYSIS  FOR 
NON  - ISOENERGETIC  SEPARATED  FLOWS 


A4.1 


APPROACH  AND  ASSUMPTIONS 


In  order  to  estabiiish  a criterion  for  reattachment  of  the  separated 
shear  layer,  a defining  set  of  conditions  must  be  imposed  upon  continuity, 
momentum,  and  energy  considerations  in  the  shear  layer  itself  and  at  the 
reattachment  point.  The  separation  bubble  behind  a backward  facing  step  has 
been  analyzed  by  Xorst  and  Page  . Their  approach  for  the  base  pres- 

sure problem  has  been  extended  in  the  present  case  to  consider  the  reattachment 
of  a separated  shear  layer  produced  by  a forward  facing  compression  surface. 

The  essential  ingredients  to  such  an  analysis  are: 

• Continuity  considerations 

• Momentum  considerations 


• Energy  considerations 

• A turning  criterion  at  .eattachment  based  upon  oblique  shock 
theory. 

The  mixing  theory  contains  the  following  assumptions: 

» The  mixing  region  is  at  a state  of  constant  pressure 

• The  length  of  the  mixing  zone  is  sufficiently  long  so  that  the 
velocity  profiles  in  the  mixing  layer  have  become  similar 


The  fluid  velocities  in  the  separation  bubble  are  small  compared 
to  velocities  in  the  mixing  zone 


• The  boundary  layer  approaching  the  point  of  separation  is  thin 


The  fluid  is  an  ideal  gas  (i.e.,  it  obeys  Boyle's  and  Gay-Lussac's 
laws) . 
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A4.2  MODEL  FORMULATION 

The  geometric  aspects  of  the  flow  separation  model  are  depicted  in 
Figure  A4.1.  Here  the  boundary  layer  approaching  the  separation  point,  a,  is 
considered  either  thin  or  nonexistant,  as  in  the  case  of  a leading  edge 
separation.  After  separation  occurs  by  whatever  mechanism  has  produced  it, 
a free-shear  layer  will  develop  between  the  external  inviscid  flow  in  Region  2 
and  the  low  velocity  flow  in  Region  4.  The  mixing  process  can  be  considered  to 
occur  at  constant  pressure  by  virtue  of  a quasi-boundary  layer  assumption  of 
zero  normal  and  axial  pressure  gradient.  As  shown,  the  abrupt  turn  of  an 
approaching  supersonic  flow  caused  by  the  separation  induces  a separation  shock 
to  form  emanating  from  point,  a,  and  will  cause  the  flow  to  turn  an  amount  612< 

At  some  distance  downstream  of  the  separation  point  (on  the  order  of 
several  approach  boundary  layer  thicknesses) , the  velocity  profiles  in  the 
shear  layer  become  similar.  The  point  of  reattachment  is  determined  by  a so- 
called  escape  criterion.  That  is,  the  point  of  reattachment  acts  as  essentially 
a stagnation  point  below  which  the  flow  in  the  shear  layer  must  return  to  the 
separation  bubble,  Region  4.  In  order  for  a steady  separation  to  exist,  the 
amount  of  flow  returned  to  the  separation  bubble  must  be  exactly  equal  to  the 
amount  cf  fluid  entrained  from  the  bubble  by  the  shear  layer.  This  condition 
also  establishes  the  existence  of  the  dividing  streamline  within  the  shear 
layer.  That  is,  all  flow  below  the  dividing  streamline  is  trapped  within  the 
separation  bubble  and,  conversely,  all  flow  outside  of  the  dividing  streamline 
never  enters  the  separation  bubble.  It  is  assumed  here,  as  is  the  usual  case, 
that  the  dividing  streamline  is  straight  between  points  a and  b. 

Upon  reattachment  of  the  dividing  streamline  to  the  ramp  surface,  the 
flow  must  undergo  a second  abrupt  change  in  flow  direction,  ^23*  This  turning 
of  the  flow  in  Region  2 at  Mach  number,  M2,  produces  a recompression  shock  as 
shown.  Hence,  the  total  amount  of  turning  produced  by  the  separation  and  re- 
attachment  shocks  must  equal  the  final  flow  direction  , 6.  Hence 


6 " 612  + 623 


(A4-1) 


A4.3  INTRINSIC  AND  REFERENCE  COORDINATE  SYSTEMS 

Because  the  mixing  equations  are  solved  subject  to  a given  simple  set 
of  boundary  conditions  independent  of  geometric  considerations  at  the  point  of 
reattachroent  and  that  determination  of  dividing  streamline  position  is  based 
upon  integral  considerations  of  continuity  and  momentum,  there  exists  the  need 
to  establish  two  sets  of  coordinate  systems.  These  two  coordinate  systems  will 
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be  termed  the  X,  Y,  or  reference  system  and  the  x,  y,  or  intrinsic  system. 

The  exact  reasons  for  establishing  two  systems  and  relating  them  will  become 
apparent  in  the  paragraphs  which  follow. 

Figure  A4.2  points  out  some  of  the  considerations  which  establish  the 
need  for  a two  coordinate  system  approach.  Figure  A4.2a  defines  the  intrinsic 
system,  x,  y,  from  which  the  velocity  profile  is  defined  based  upon  the  edge 
conditions  in  regions  2 and  4.  The  profile  independent  variable  is  defined  as, 
n = oy/x,  in  the  intrinsic  system  where  a is  the  turbulent  mixing  parameter. 

The  extent  of  jet  spreading  is  defined  by  lines  a and  e in  Figure  A4-2b. 

Figure  A4.2b  shows  the  relation  of  the  reference  system  and  the  intrinsic  system. 
The  reference  system  is  orthogonal  and  aligned  with  the  u^  direction.  Introduc- 
tion of  the  profile  defined  in  the  intrinsic  system  shown  as  line  d is  required 
since  the  profile  is  described  in  terms  of  the  variable,  n.  It  is  necessary, 
therefore,  to  define  the  displacement  of  the  Y and  y coordinates,  namely,  y^. 

It  will  be  shown  that  this  displacement  between  the  direction  (line  b,  refer- 
ence system)  and  the  intrinsic  direction  (line  d)  can  be  determined  from  the  in- 
tegral relations  defining  the  location  of  the  dividing  streamline,  c.  Hence,  the 
two  coordinate  systems  coincide  at  x = X = 0 and  elsewhere  are  related  by 


X s x 


(A4-2a) 
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A4.4  VELOCITY  AND  TOTAL  TEMPERATURE  PROFILES 

The  turbulent  equations  of  motion  and  energy  in  two-dimensional  space  are 
given  by. 


pu 


(A4-3) 
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(A4-4 ) 


Equation  (A4-4)  is  derived  based  on  the  assumption  of  unity  Prandtl  number.  Hence, 
it  can  be  seen  that  Equations  (A4-3)  and  (A4-4)  are  analogous  if  the  boundary  condi- 
tions are  made  compatible.  Equation  (A4-3)  can  be  further  simplified  by  a linear- 
ization scheme  described  in  Reference  A17  and  reduced  to  the  following  form. 
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A transformed  set  of  coordinates  is  introduced  as  described  in  Reference-  15A 
to  yield  an  equation  similar  to  the  heat  conduction  equation.  The  final  solu 
tion  to  this  equation  is  given  by 


♦ 


“ l(l  + erf  n) 


(A4-6) 


where 


(A4-7 ) 


o = turbulent  mixing  parameter  = 12 
(typical  turbulent  value) 


and 


erf (n) 


(A4-8) 


The  energy  equation  may  be  treated  in  a similar  manner  and,  since  Equations  (A4-3) 
and  (A4-4)  are  analogous,  the  result  is 


(A4-9) 


or 
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The  static  temperature  ratio  will  be  given  by 
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where  the  Crocco  number  in  Region  2 is  defined  by 


C2  = u 


max  2/  (y  -1)  + M2 


(A4-14) 


and  temperature  ratios  are  defined  by  isoenergetic  expansion  from  the  local  stag- 
nation conditions. 

A4.5  INTEGRAL  MOMFNTUM  AND  CONTINUITY  CONSIDERATIONS 

Figure  A4-3  illustrates  the  two  control  volumes  of  interest,  CV1  and 


A4 . 5 . 1 Momentum  Considerations 

Since  the  mixing  process  occurs  at  constant  pressure  and  since  no  other 
extcnral  forces  act  on  CV1  in  Figure  A4-3a,  conservation  of  momentum  is  defined 
by  the  fact  that  the  stream  momentum  remains  cor  tant. 
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Neglecting  any  momentum  deficiency  in  the  thin  approaching  boundary  layer  and 
noting  that 


at  (1) 
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n = 21  = 2l 

1 x x 
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at  (2) 


Y = y - y. 
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dY  = dy 
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Then,  at  (1) 


and  at  (2) 


Y Y y 

/R  /r  >R 

/ pu 2 dY  - p2u|  I dY  = p2u|  / dy 
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pu2dY  = / pu2dy  = / pu2dy  + p2u2Jyn 


Hence 


/"pu2dy  = p2u22  A = P2u22ym 


i u?2yP. = y pu*d* + p2u22yn 


a. 

r,  = I — <>2dy  + r, 

R J P2  * «** 


Now,  also  for  constant  pressure  mixing  of  an  ideal  gas 


T i - r - 

_P _ = _2  _ "2 

o,  T x - c 

2 A - *_2  > 


so  that 


where 


= nR  " W 


I ( ,)  = (!-  C 2)  f 

2 *•  2 J X - C 24? 


Hence,  the  distance  between  the  intrinsic  and  reference  coordinate  systems  is 
given  by  Equation  (A4-26) . 
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A4.5.2  Continuity  Considerations 

The  escape  criterion  described  previously  to  define  the  dividing  stream- 
line is  used  to  determine  the  location  of  the  dividing  streamline  with  respect 
to  the  intrinsic  coordinate  system.  Figure  A4-3b  illustrates  the  control  vol- 
ujije t CV2 , above  the  dividing  streamline.  Hence, 
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In  dimensionless  form  eliminating  the  density  ratio  and  neglecting  the  massflow 
deficiency  in  the  approaching  boundary  layer,  we  have 
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Substituting  Equation  (A4-26)  and  noting  thac 


n.  1 - C 
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yields 


1 1 (n j ) = ^l^R^  “ I2^nR^ 
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where 
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Therefore,  if  1^  and  I,  are  evaluated  up  to  the  limit  of  the  inviscid  boundary 
(say  = 3)  then  the  position  of  the  dividing  streamline  in  terms  of  can  be 
deduced  from  the  value  of  I j (n j ) - 

A4.6  MASS  ADDITION 

If  mass  is  added  or  subtracted  from  the  separate-*  zone  as  depicted  in 
Figure  A4-4  in  order  to  retain  continuity  an  equal  amount  of  fluid  would  have 
to  leave  or  be  added  tc  the  separated  zone  if  the  separation  were  to  remain  stable. 
That  is,  the  fluid  being  added  would  ultimately  have  to  escape  through  the 
plane  of  reattachment.  This  would,  therefore,  cause  the  point  of  reattachroent 
to  shift  from  the  j streamline  to  a new  streamline  defined  as  the  S streamline. 
This  then  becomes  simply  a displacement  of  the  attaching  streamline  by  an  amount 
which  allows  the  injected  flow  to  escape. 


W - w - Gib 


(A4-36 ) 


where 


Gib  " 


(A4-37) 


Here  is  the  dimensionless  blowing  (G^  > 0)  or  bleed  (G^  < 0)  rate  from  the 
separated  zone,  is  the  corresponding  massflow  per  unit  width,  and  l is  the 
length  of  the  separation. 

Inspection  of  Equation  (A4-36)  shows  that  if  flow  is  added  or  subtracted 
from  the  separated  zone  the  velocity  at  reattachment  will  be  smaller  or  greater 
respectively,  than  the  j streamline  velocity.  The  following  paragraph  will 
detail  how  this  will  affect  the  allowable  turning  at  reattachment.  It  will  be 
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shown  in  Section  A6.2  how  in  the  case  of  the  reattachmant  condition  being 
violated, it  may  be  possible  to  relate  the  rate  of  bubble  growth  to  the  para- 
meter, G^fa.  In  this  instance  the  parameter  is  to  be  thought  of  as  the 
driving  potential  for  bubble  growth. 

REATTACHMENT  CONDITIONS 

The  Mach  number  on  the  reattaching  streamline  would  be  given  by 
: for  Gib  = 0 


u . 

M.  U 

M — 1 = M — 

2 m2  2 fr 


r-j 

= M,  -1 
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: and  for  G^b  ? 0 
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MS=M2  = 
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In  order  to  determine  the  conlitions  after  reattachment  the  following  approach 
was  taken.  If  the  Mach  number  on  the  reattaching  streamline  is  subsonic  then 
just  before  reattachment 


fc2r  _ /.  Y _ 1 ..  j 


("r  < ') 


(A4-38) 


where  r = j or  s 

and  since  the  reattachment  point  is  considered  essentially  a stagnation  point 


P,  = P. 
3r  t. 


(at  reattachment) . 
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If,  on  the  other  hand,  the  reattachment  Mach  number  is  greater  than  1.0  then 
the  recompression  shock  which  forms  within  the  shear  layer  must  bend  nearly 
normal  to  the  approaching  flow.  Then  the  stagnation  pressure  behind  the  recom- 
pression shock  and  at  the  stagnation  point  will  he  given  by  the  Rayleigh  pitot 
formula. 
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P,  * P (at  reattachment) 
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In  both  instances  the  static  pressure  level  behind  the  reattachment  shock  is 

assumed  to  be  equal  tc  the  reattachment  pressure.  This  technique  was  used  with 

. „ (A15) 

success  by  Korst 

The  allowable  turning  at  reattachment  is  then  determined  by  the  oblique 
shock  relations  at  reccxnpression.  The  recompression  shock  angle,  e23' 
respect  to  the  flow  direction  in  region  2 is  given  by. 


023  = Sin 


U-  + 


y + i 


(A4-42) 


where 


: 

p23  yM22 


(A4-43) 


and  the  flow  is  turned  an  amount,  ^23'  given  by 


5 = cot  /tan  G 

23  1 M22sin29  - 


JZj  M2?  - 1 ( . 

sin20  - 1 ^ 


The  flow  on  the  reattachment  streaml.  .e  will  turn  an  amount 


(A4-44) 


6r  = 623  + Oj.  •••  (r  = j or  s) 


(A4-45) 


where 


ar  = ten 


■(^) 


(A4-46) 
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SECTION  A5 

FLOW  SEPARATION  INTERACTIONS 
A5.1  FREE  AND  FORCED  INTERACTIONS 

Two  types  of  turbulent  boundary  layer  separations  may  be  defined  as  those 
which  are  free  interactions  and  those  which  are  forced  interactions.  A free 
interaction  may  be  defined  as  one  in  which  the  separation  geometry  is  not 
dependent  upon  any  downstream  effects  or  what  produced  the  separation.  Figure 
A5-la  illustrates  some  of  the  salient  features  of  the  free  interaction.  The 
flow  approaching  the  point  of  separation  has  a well  defined  turbulent  boundary 
layer.  The  flow  deflection,  ^12'  at  separation  and  the  rise  to  plateau  pressure 
level  are  defined  by  a delicate  balance  of  pressure  and  momentum  forces  at  the 
point  of  separation.  Therefore,  the  cause  of  the  separation,  whether  it  be  a 
ramp,  curved  surface  or  forward  facing  step,  is  of  no  consequence  in  determining 
the  angular  geometry  at  separation. 

Suppose,  however,  the  distance  from  the  corner  to  the  leading  edge  is 
significantly  less  than  the  typical  length  associated  with  a f ree- interaction , 
or  suppose  than  the  separation  is  forced  to  occur  at  or  near  the  leading  edge, 
as  shown  in  Figure  A5-lb.  It  will  no  longer  be  true  that  the  geometry  of  separa- 
tion, 5j2'  be  independent  of  downstream  influences.  In  fact  the  angle, 

c12,  may  be  totally  dependent  upon  the  downstream  reattachment  condition. 

A5.2  FREE  INTERACTION  MODEL 
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In  the  instance  of  a free  interaction  separation  the  geometry  at  separa- 
tion was  deduced  from  the  empirical  plateau  pressure  rise  coefficient  found  in 
Reference  Alj. 
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12 


1 + 10  M. 


(A5-1) 


This  correlation  reflects  the  fact  that  for  fully  turbulent  separations  the 
plateau  pressure  rise  coefficient  is  largely  independent  of  Reynolds  number. 
Th«i  separation  angle  with  respect  to  the  approaching  flow  was  determined  using 
the  locally  two-dimensicnal  oblique  shock  relations  given  by  equations  similar 
to  A4-42  and  A4-44.  The  Mach  number  in  region  2 is  given  by 
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Hence,  the  free  interaction  plateau  pressure  coefficient  is  sufficient  to  specify 
all  required  conditions  to  proceed  with  the  reattachment  analysis.  The  total 
amount  of  turning  sustainable  by  an  asymptotic  free  interaction  is  given  by 

0 = 6i2  + 623  (A5-3) 

where  6 is  obtained  from  Figure  A5-2  and  623  is  obtained  by  th«i  methods  detailed 
in  Section  A4. 


AS. 3 FORCED  INTERACTION  MODEL 

The  forced  interaction  model  is  similar  in  concept  to  a model  applied 
to  forced  laminar  separations  occurring  right  at  the  leading  edge  found  in 
Reference  A19.  In  the  present  instance  the  separation  is  not  forct  j to  the 
leading  edge  but  instead  to  the  vicinity  of  the  tangent  point  between  the  sphere- 
like nose  and  the  first  conical  surface.  The  calculation  proceeds  as  follows: 

• At  a fixed  freestream  Mach  number,  Mq  = , approaching  the  point  of 

separation  the  flow  is  allowed  to  sustain  some  specified  amount  of 
initial  turning,  5^,  caused  by  the  forced  separation. 

• The  conditions  in  region  2 are  obtained  by  either  two-dimensional 
oblique  shock  theory  or  conical  shock  theory  ignoring  the  entropy 
gradient  due  to  the  rounded  nose. 

• The  reattachment  conditions  are  determined  by  the  methods  outlined 
in  Section  A4  to  obtain  the  turning  at  re  attachment , <$23* 

• The  maximum  second  flare  angle  S for  reattachment  on  the  face  of 
the  second  flare,  is  then 

9max  = ^12  + 623  (A5-4) 

The  conical  shock  theory  used  is  discussed  in  Reference  A20.  The  two 
dimensional  reattachment  criteria  developed  in  Section  A4  should  be  nearly 
correct  for  most  important  ax j symmetric  cases  considered  in  the  context  of  this 
study  if  one  accepts  the  error  function  profile  as  being  nearly  correct.  The 
reason  for  this  is  that  at  the  reattachment  point  for  large  values  of  S there  is 
only  a small  percentage  radial  change  from  the  innermost  to  outermose  edges  of 
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the  shear  layer.  Hence,  a radial  integration  of  the  pertinent  integrals 
will  not  be  significantly  different  than  the  two  dimensional  integrals.  In  the 
case  of  axisymmetric  base  flows  the  reattachment  is  close  to  the  centerline  and 
the  radial  effect  may  be  important  as  discussed  in  Reference  A21. 

Note  that  for  the  configurations  considered  in  thi * study  there  is  no 
effect  of  the  physical  surface  angle  9^  in  the  forced  interaction  model.  In 
terms  of  the  free  interaction  model  if  the  separation  occurs  downstream  of  the 
leading  edge  then  9^  enters  the  calculation  by  virtue  of  relating  the  freestream 
Mach  number,  Mq,  to  the  Mach  number  ahead  of  the  point  of  separation  as  shown 
in  Figure  A5-3a.  As  will  be  discussed  in  Section  A6  it  is  concluded  that  for 
the  configurations  considered  in  this  study  that  all  separations  are  of  the 
forced  type.  The  questions  posed  serve  to  define  the  following  procedures: 

• Determine  for  the  particular  configuration  of  interest  S) 

whether  an  asymptotic  attached  solution  is  possible. 

• If  the  free-interaction  solution  is  not  possible  then  determine  if 
a forced  interaction  solution  is  possible. 

• If  neither  is  possible  then  the  separation  may  try  either  to  attach 
near  the  shoulder  or  force  the  flow  to  become  unstable. 
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SECTION  A6 
RESULTS 

A6.1  EFFECT  OF  NON-ISOENERGETIC  FLOW 

The  relations  described  in  Sections  A4  and  A5  were  coded  in  Fortran  IV 
and  programmed  on  an  UNIVAC  1108  digital  computer.  In  order  to  determine  if  the 
effects  of  heat  transfer  to  the  separated  region  represented  a significant 
effect  upon  the  asymptotic  reattachment  model  a set  of  sample  calculations  was 
performed.  Figure  A6-1  presents  the  results  of  this  investigation  for  a 2-D 
free-interaction.  Clearly,  relatively  large  differences  in  total  temperature 
do  not  alter  the  amount  of  turning  at  reattachment,  significantly.  Therefore, 
the  effects  of  non-iscenergetic  flow  was  not  considered  further. 


it 


£ 


A6.2  EFFECT  OF  BLEED  FLOW 

Figure  A6-2  represents  a study  of  the  effect  o'  bleeding  flow  from  the 
separated  region  of  a 2-D  free-interaction.  Clearly,  a significant  increase  in 
the  total  amount  of  turning  can  be  realized  if  flow  is  removed  from  the  separa- 
tion zone.  Since,  in  the  present  application  no  flow  removal  can  be  realized 
one  may  ask  what  significance  can  be  placed  upon  this  study.  The  following 
interpretation  is  offered.  If  the  separation  cannot  remain  stable  as  either  a 
free  interaction  or  as  a forced  interaction  an  unsteady  growth  of  the  bubble  will 
occur.  The  collapsed  bubble  will  begin  to  grow  from  the  corner.  As  an  example 
of  how  Gib  may  be  used  to  estimate  the  rate  of  bubble  growth  consider  the 
following.  Say  the  approach  Mach  number  at  the  beginning  of  bubble  growth  is 
= 4.0,  and  say  8 - 42°.  Examination  of  Figure  A6-2  shows  that  if  G^b  =-.20, 
that  is,  if  that  quantity  of  flow  is  removed  the  separation  will  be  steady  and 
stable.  However,  if  the  outflow  valve  is  shut  off  so  to  speak  then  instantaneously 
the  bubble  will  accumulate  that  amount  of  flow  which  is  not  removed  by  bleed. 

Hence,  the  rate  of  a non-dimensional  bubble  growth  is  G^b  = .20.  If  on  the 
ether  hand  8 exceeds  the  allowable  limiting  value  at  = 4.0  namely  6max  = 47.5° 
then  the  rate  of  bubble  growth  is  interpreted  as  G..  = .36. 


A6.3  FREE  AND  FORCED  INTERACTION  RESULTS 

Two  sets  of  forced  interaction  results  were  generated  and  are  shown  in 
Figures  A6-3  and  A6-4.  The  results  shown  in  Figure  A6-3  are  for  a two-dimensional 
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A6-4.  Forced  Leading  Edge  Interaction  Result 

• Conical  Separation 
t 2-0  Shear  Flow 

• Y = 1-4 
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forced  interaction.  The  curved  labeled  610  • „ represents  the  free- interaction 

12  mm 

solution.  The  results  shown  in  Figure  A6-4  are  for  a forced  separation  initiated 
with  a conical  shock.  The  results  of  Figures  A6-3  and  A6-4  do  not  differ 
greatly  although  the  allowable  amount  of  turning  is  somewhat  greater  for  the 
axisymmetric  case.  As  mentioned  previously  the  reattachment  results  for  Figure 
A6-4  were  performed  using  the  two-dimensional  model;  however,  the  radial  effect 
at  reattachment  should  be  small  for  the  larger  values  of  8 


SECTION  A7 


CONCLUSIONS 

A supersonic  asymptotic  separation  and  reattachment  model  has  been  formu- 
lated to  predict  whether  the  flow  separation  on  a biconic  surface  will  be  stable 
or  unstable.  The  model  is  formulated  on  the  basis  of  the  initial  separation 
being  either  a free  interaction  or  a forced  interaction.  The  reattachment  is 
based  on  a free-shear  flow  analysis.  If  the  separation  is  a forced  interaction 
the  criteria  for  maximum  allowable  turning  is  independent  of  the  first  cone 
angle.  The  effects  of  heat  transfer  to  the  separated  zone  are  small  in  relation 
to  the  allowable  turning  angle.  The  effect  of  bleed  flow  has  been  related  to 
the  growth  rate  of  the  separation  bubble  for  the  unstable  cases.  The  forced 
separation  cases  illustrate  that  the  maximum  amount  of  turning  can  be  signif- 
icantly greater  than  for  the  f ree-interacticn.  The  results  of  the  present 
study  indicate  that  the  maximum  allowable  turning  for  a conical  reattachnent 
is  about  55°  at  Mach  5. 

It  seems  clear  that  the  weakness  of  the  rodel  lies  in  its  inability  to 
predict  occurences  after  the  sonic  limit  has  been  reached.  Once  the  flow 
becomes  locally  setae  a closer  coupling  to  the  inviscid  flow  would  be  necessary 
to  predict  the  reattachment  criteria.  This  is  clearly  beyond  the  scope  of 
the  model  presented  herein. 


A7-1 


mmrn&Siz atm 


i k-'Ka-'  v ■ 


REFERENCES 


FOR  APPENDIX  A 


Wood,  C.  J.,  "Hypersonic  Flow  over  Spiked  Cones"  J.  of  Fluid  Mechanics, 

Vol.  12,  Part  4 April  1962. 

Maull,  D.  J. , "Hypersonic  Flow  over  Axially  Symmetric  Spiked  Bodies", 

J.  of  Fluid  Mechanics,  Vol.  8,  Part  4.  Aug.  1960. 

Mair,  W.  A.  Phil.  Mag.  1952,  43,  695. 

Zukoski,  E.  E.,  "Turbulent  Boundary-Layer  Separation  in  Front  of  a Forward* 
Facing  Step"  AIAA  Journal,  Vol.  5,  No.  10,  Oct.  1967. 

Love,  E.  S.,  "Pressure  Rise  Associated  with  Shock-Induced  Boundary-Layer 
Separation",  NACA  TN  3601,  Dec  1955. 

Settles,  G.  S.,  and  Bogdor.off,  S.  M. , "Separation  of  a Supersonic 
Turbulent  Boundary  Layer  at  Moderate  to  High  Reynolds  Numbers."  AIAA 
Paper  73-666,  July  1973. 

Kuehn,  D.  M. , "Experimental  Investigation  of  the  Pressure  Rise  Required 
for  the  Incipient  Separation  of  Turbulent  Boundary  Layers  in  Two 
Dimensional  Supersonic  Flow”  NASA  Memo  1-21-59 A,  Feb.  1959. 

Law,  C.  H.,  "Supersonic,  Turbulent  Boundary-Layer  Separation  Measurements 
at  Reynolds  Numbers  cf  107  to  10'’",  AIAA  Paper  73-665,  July  1973. 

Roshko,  A.,  and  Thomke,  G.  J. , "Supersonic,  Turbulent  Boundary-Layer 
Interaction  with  a Compression  Corner  at  Very  High  Reynolds  Number". 

Chapman,  D.  R. , Kuehn,  D.  M. , and  Larson,  H.  K.,  "Investigation  of 
Separated  Flows  in  Supersonic  and  Subsonic  Streams  with  Emphasis  on 
the  Effect  of  Transition",  NACA  TN  3863,  Mar.  1957. 

Sterrett,  J.  R, , and  Emery,  J.  C.,  "Experimental  Separation  Studies  for 
Two-Dimensional  Wedges  and  Curved  Surfaces  at  Mach  Cumbers  of  4.8  to 
6.2,"  NASA  TND-10M  Feb.  1962. 

TcJisco,  A.,  and  Reeves,  B.  L. , "Turbulent  Boundary  Layer  Separation  and 
Reattachment  at  Supersonic  Speeds". 

Spaid,  F.  W. , "Cooled  Supersonic  Turbulent  Boundary  Layer  Separated  by 
a Forward  Facing  Step,"  AIAA  Journal,  Vol.  IV,  No.  8.,  Aug.  1972. 

Kuehn,  D.  M. , "Turbulent  Boundary-Layer  Separation  Induced  by  Flares  on 
Cylinders  at  Zero  Angle  of  Attack,"  NASA  TRR-117. 

Korst,  H.  H. , "A.  Theory  for  Base  Pressures  in  Transonic  and  Supersonic 
Flow,"  J.  o£  Applied  Mechanics,  Dec.  1956. 

Page,  R.  H.,  and  Korst,  H.  H.,  "Non-Isoenergetic  Compressible  Jet 
Mixina  with  Consideration  of  its  Influence  on  the  Base  Pressure  Problem," 
presented  at  the  Fourth  Midwestern  Conference  on  Fluid  Mechanics, 

Purdue  1955. 


Al7. 


Korst,  H.  H. , Page,  R.  H.,  and  Childs,  M.  E.,  "A  Theory  for  Base  Pressures 
in  Transonic  and  Supersonic  Flow,"  ME  TN  392-2,  OSR-TN-55-39 , U of 
Illinois,  March  1955. 

A18.  D.hm,  T.  J.,  "A  Comprehensive  Analytical  Procedure  for  the  P3.*';ormance 
Prediction  of  Rocket  Thrust  Vector  Control  with  Gaseous  Secondary 
Injection,"  Aerotherm  Report  67-11,  1967. 

A19.  Brower,  M.  3.,  "Leading-Edge  Separation  of  Laminar  Boundary  Layers  ir. 

Supersonic  F.ow",  J.  of  the  Aerospace  Sicences , Dec.  1961. 

A20.  Moretti,  G.  "Inviscid  Flow  Field  Past  a Pointed  Cone  at  an  Angle  of 
AttacL  Pvrt  I"  GASL  Tech.  Report  577,  Dec.  1965. 

A21.  Mueller,  T.  J.,  “Determination  of  the  Turbulent  Base  Pressure  in  Super-  * 

sonic  Axisymmetric  Flow,"  J.  Spacecraft,  Jan.  1968.  I 


